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An  orthotropic  right  cylindrical  shell  is  analyzed  when  subjected  to  a  dis- 
continuous, finite  length  pressure  load  moving  in  the  axial  direction  at  constant 
velocity.    The  analysis  utilizes  linear,  small  deflection  shell  theory  which  in- 
cludes the  effect  of  axial  and  circumferential  prestress,  transverse  shear  de- 
formation, and  external  radial  damping. 

The  problem  is  solved  using  Fourier  transforms,  and  the  inverse  Fourier 
integrals  are  evaluated  for  the  radial  deflection,  axial  deflection  and  rotation 
by  expanding  the  Characteristic  Equation  in  partial  fractions  and  using  complex 
contour  integration.    By  studying  the  discriminant  of  the  undamped  characteris- 
tic equation  the  load  velocities  which  give  repeated  roots  are  determined.    The 
loci  of  these  load  velocities  separate  regions  in  which  the  form  of  the  displace- 
ment solutions  differ.    The  behavior  of  these  load  velocity  loci  is  studied  for 
variations  in  the  three  nondimensionalized  material  moduli ,  the  thickness-to- 
radius  ratio,  the  axial  prestress,  and  the  circumferential  prestress. 
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By  tracing  the  root  loci  of  the  undamped  characteristic  equation  and  by  in- 
spection of  the  displacement  expressions,  it  is  determined  that  there  are  five 
critical  load  velocities  (velocities  at  which  the  displacement  becomes  unbounded) 
for  the  specific  example  of  an  isotropic  shell.    An  increase  of  the  load  velocity 
above  the  bar  wave  speed  produces  a  deflection  mode  which  is  predominantly 
axial. 

The  deflection  response  is  investigated  for  numerous  combinations  of  load 
speed,  material  properties,  length  of  pressure  load,  axial  and  circumferential 
prestress,  and  radial  damping.    The  axial  prestress  has  a  significant  effect  on 
the  first  critical  velocity  of  the  cylinder;  initial  compression  tends  to  lower 
the  velocity .    Circumferential  prestress  has  no  pronounced  effect  on  the  critical 
load  speeds  but  does  influence  the  response  at  higher  velocities.    Variation  of 
material  properties  was  found  to  cause  a  rapid  change  in  deflection  response. 

Through  superposition,  the  variation  of  pressure  load  length  can  be  utilized 
to  approximate  the  response  to  any  desired  pressure  load.    Examples  of  this 
application  are  demonstrated.    A  comparison  of  stresses  and  deflections  against 
those  predicted  by  the  Timoshenko  thin  shell  theory  is  shown  for  a  static  load . 

All  of  the  above  numerical  work  was  done  using  dimensionless  parameters 
which  can  be  applied  to  thin  shells  in  general.    The  calculations  were  done  utiliz- 
ing a  computer  program  developed  from  this  research  for  the  calculation  of  de- 
flections and  stresses  in  the  shells.    The  program  is  written  in  Fortran  and  is 
operable  on  the  General  Electric  Company  Mark  II  time  sharing  service . 
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CHAPTER  I 
INTRODUCTION 

Statement  of  the  Problem 

One  of  the  most  commonly  used  geometries  for  structural  application  is  the 
right  circular  cylindrical  shell.    This  is  particularly  true  in  the  aerospace  field 
and  in  undersea  exploration  vehicles.    In  many  aerospace  applications  the  cylin- 
drical shell  serves  as  the  primary  load  carrying  member  for  the  rocket  system 
and  performs  simultaneously  as  a  portion  of  the  pressurized  fuel  tank.    In  under- 
sea applications  the  quest  for  greater  depth  range  has  brought  about  many  re- 
finements in  structural  optimization  techniques.    A  result  of  the  many  stringent 
requirements  being  placed  upon  structural  systems  has  resulted  in  two  areas  of 
rapid  advancement:    new  material  technologies  and  more  sophisticated  analysis 
techniques . 

The  material  technologies  for  advanced  design  applications  have  in  many 
cases  moved  away  from  the  isotropic  materials  and  are  utilizing  orthotropic  and 
anisotropic  materials  to  satisfy  the  demanding  requirements  for  more  efficient, 
lighter  weight  vehicles.    Studies  such  as  that  reported  in  Reference  (1)*  have 
shown  that  there  is  indeed  an  incentive  for  the  application  of  these  advanced 
technologies. 

Until  recent  years  the  mathematical  complexity  encountered  when  approach- 
ing the  dynamic  analysis  of  shells  has  been  so  formidable  that  few  results  were 
available  for  design  applications . 


♦Denotes  entries  in  the  Bibliography, 


The  new  technology  demands  mentioned  previously  have  brought  a  response 
from  the  analysts  in  the  past  five  to  ten  years  and  some  of  the  more  idealized 
dynamic  shell  problems  have  been  investigated.    The  problem  of  particular  in- 
terest in  this  work  is  that  of  a  thin  orthotropic  cylindrical  shell  subjected  to  an 
axisymmetric  pressure  load  moving  in  the  axial  direction.    It  is  necessary  to 
consider  refinements  to  the  theory  such  as  the  transverse  shear  deformation, 
axial  inertia,  and  rotatory  inertia  effect  so  that  the  higher  load  velocities  may 
be  investigated. 

It  is  known  that  axial  prestress  has  an  influence  on  results  in  dynamic 
analyses.    In  this  work  the  effect  of  axial  as  well  as  circumferential  prestress 
is  investigated.    The  specific  loading  considered  will  be  a  constant  pressure 
pulse  finite  in  both  magnitude  and  distance  which  moves  along  the  cylinder 
at  velocity  V.    The  shell  theory  utilized  is  linear,  assuming  small  deflections, 
and  by  superposition  it  is  possible  to  investigate  the  effect  of  various  pressure 
pulse  shapes.    External  radial  damping  is  also  included. 

Specific  Goals  of  This  Research 

The  major  goal  of  this  research  was  to  obtain  a  solution  for  the  deflections 
and  stresses  associated  with  the  problem  outlined  above.    The  secondary  goal, 
although  perhaps  not  secondary  in  importance  to  those  interested  in  utilizing  the 
results,  was  that  of  developing  a  computer  program  for  the  calculation  of  deflec- 
tions and  stresses  in  the  cylinder.    Finally,  the  calculation  and  presentation  of 
the  effect  of  the  many  parameters  included  in  the  analysis  conclude  the  goals  to 
be  reached  in  the  study . 


Review  of  Previous  Work 

A  review  of  the  early  work  on  the  response  of  a  cylindrical  shell  to  a  moving 
load  is  given  by  Jones  and  Bhuta  (2).     Until  the  work  by  Nachbar  (3),  who  con- 
sidered the  dynamic  response  of  an  infinitely  long  cylindrical  shell  to  a  semi- 
infinite  step  pressure  load,  the  axial  inertia  was  not  considered.    Nachbar 
included  the  axial  inertia  effect  and  also  assumed  an  external  damping  effect. 
However,  due  to  the  damping  included,  the  first  resonance  condition  was  missed, 
Jones  and  Bhuta  solved  the  problem  with  a  ring  load  moving  on  an  infinitely  long 
cylinder  but  did  not  include  the  transverse  shear  effect. 

Other  contributions  were  made  by  Reismann  (4)  who  included  the  effect  of 
axial  prestress,   which  was  significant  as  had  been  found  in  his  work  on  plate 
strips  (5).    Hegemier  (6)  studied  the  stability  problem  for  a  large  class  of  con- 
stant velocity  moving  loads  but  limited  the  velocity  range  to  that  lower  than  the 
first  critical. 

All  of  the  work  previewed  above  was  done  for  an  isotropic  material.    More 
recently  Herrmann  and  Baker  (7)  solved  the  problem  of  a  moving  ring  load  on  a 
cylindrical  sandwich  shell  of  infinite  length.    Numerical  results  were  presented 
for  a  core  material  which  is  assumed  to  have  material  damping.    Also,  the 
problem  of  a  ring  load  moving  on  a  viscoelastic  cylinder  was  solved  by  Tang  (8) . 

Contributions  of  This  Work 

The  following  contributions  are  believed  to  be  original  with  this  work. 

1 .  Analysis  of  orthotropic  monocoque  cylindrical  shells  including  trans- 
verse shear  deformation,  axial  and  rotatory  inertia,  radial  damping, 
circumferential  and  axial  prestress  under  a  finite  length  step  load. 

2.  Presentation  of  the  forms  of  the  solution  with  no  damping  for  the  seven- 
dimensional  space  whose  coordinates  are  the  thickness-to-radius  ratio, 


the  three  material  property  ratios,  axial  prestress,  circumferential 

prestress,  and  the  load  velocity  parameter. 

Results  indicating  the  effect  of  a  finite  length  pressure  pulse,  and  the 

capability  to  approximate  any  load  shape  through  superposition. 

Indication  of  the  effect  of  prestress  on  the  critical  velocities  of  an 

orthotropic  monocoque  shell. 

Results  which  show  the  effects  of  external  damping  throughout  the  load 

velocity  range. 


CHAPTER  n 
GOVERNING  EQUATIONS  OF  MOTION 

General  Equations 

A  cylindrical  shell  of  thickness  h  and  mean  radius  R  is  referred  to  the  co- 
ordinate system  shown  in  Figure  2.1. 


Figure  2.1.    Cylindrical  Coordinate  System 

Coordinate  x  is  measured  along  the  shell  axis,   e  along  the  circumference 
and  z  is  perpendicular  to  the  middle  surface.    The  unit  vectors  tangent  to  the 
coordinate  lines  at  a  point  (x,  e,  z)  are  designated  by  K  ,  K      K  .    The  dis- 
placements in  these  three  directions  are  u   ,  u      and  u    respectively.    It  is 


assumed  that  the  displacements  can  be  represented  by  the  linear  relationship  in 
terms  of  z 

u^(x,  9,  z,  t)     =     u(x,  e,  t)  +  z>p{x,  e,  t) 

u^(x,  e,  z,  t)     =     v(x,  9,  t)  +  ZT7(x,  e,  t) 

u^(x,  0,  z,  t)     =     w(x,  9,t)  +  z^^(x,  9,  t)  (2-1) 

where  u,  v,  w  are  the  displacements  on  the  middle   surface  (z  =  0);  t  denotes 
time ;  ip  and  tj  are  the  rotations  of  a  line  perpendicular  to  the  normal  surface  in 
the  x-z  and  9-z  planes  respectively,    ip    is  referred  to  as  the  thickness  stretch. 

Equations  (2-1)  require  that  all  straight  lines  normal  to  the  middle  surface 
of  the  shell  before  deformation  remain  straight  after  deformation.    This  is  a 
good  approximation  if  the  shell  is  thin. 

Herrmann  and  Armenakas  (9)  derived  a  linearized  theory  for  the  motion  of 
isotropic  cylindrical  shells  subjected  to  a  general  state  of  initial  stress  by  as- 
suming the  final  state  of  stress  is  reached  by  passing  through  an  intermediate 
state,  the  state  of  initial  stress.    Subtracting  the  initial  equilibrium  equations 
from  the  non-linear  equations  of  motion  and  then  linearizing  by  disregarding  all 
non-linear  terms  involving  the  additional  stresses,  the  linearized  equations  of 
motion  for  a  shell  under  initial  stress  are  obtainable. 

Following  this  procedure,  Baker  and  Herrmann  (10)  derived  a  linearized 
set  of  equations  for  the  motion  of  orthotropic  shells.    Assuming  an  orthotropic 
cylindrical  monocoque  shell  is  under  initial  lateral  pressure  p.,  an  axial  tension 
T,  and  is  subjected  to  a  time  dependent  radial  load,  the  five  equations  of  motion 
have  the  form: 

X    ,   1         9x    ,  „  9  u    ,     N     9  u        N    9w  u  .    I    ,  /o   o\ 

IF^R-gr^T^-*-;^^  -R  gT    =    ^  V  ^  R  ^  (^-2) 


1^  

R    dd    "      Ox 


xp 


«  8x^     R^^^      a/ 


hn   V 


I  ••       N  r^o  aw  ^, 


(2-3) 


N^^aO^^^aO,         93^    _N/    ^av     afw  ^  av 


+  q'(x,  0,  t)     =     ^w  +  p  hw 


(2-4) 


aM        .   aM„ 

X  ^  _1  Ox 


ax      R      ax 


-Qx-Nf  (.^If) 


I'^'  +  ^u 


(2-5) 


R    ae        ax       ^e 


5[-('4)if^*o-^ 


(2-6) 


where 


h 
2 


h 
2 


I    = 


1,3 
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p.R     1  +  r^  for  external  pressure 


N     = 


\^['-£ 


(2-7) 


for  internal  pressure 


(2-8) 


These  equations  include  the  effect  of  external  radial  damping,  axial  and 
rotatory  inertia,  and  transverse  shear  deformation. 


8 

The  elongations,  shears,  and  rotations  have  been  assumed  small  in  com- 
parison with  unity.    The  strain-displacement  relations  are  therefore  taken  in 
the  form  of  Hooke's  law. 

^     8u 

^x         ax 

^e        R  +z  \9e      ^J 

"^xe        R  +  z    ae      ax 

au  ^  aw 

"^xz        az      ax 

9v  ^     1     /aw      \  -o  ON 

^0z  ^  ai  "-RTiy-der-  V  <2-^) 

The  stress  strain  relationships  are  assumed  in  the  following  form 

a        =     Ex    C     +  Ey     C. 
X  ^o     x  ^o      9 


=     E„     C     +  Eq     C. 
'^o      X  "o      6 


''xz  *^xzo  \z  (2-10) 


Integrating  the  stresses  through  the  shell  thickness  the  z  variable  is  eliminated 
and  the  stress  and  moment  resultants  are  obtained. 


N 


^x  ax  ^  ^  37  ^  ^  (^^  ^  90  ; 


N 


'■e    .  ''e 


R 


w   + 


R^ 


8v 

do 


_i    92.  +   E     ^ 
j^2     de  V  9x 


N 


xe 


D  P 

G      ^  +  _2ii  M  +  _Ji£  9u 
xe  ax        R    ax        R    dd 


N, 


ex 


G       9v        5ii   9i   ^   /5s£   +  "x0\  9u 

xe  ax       r2   ap      \  R        ^3  j  39 


M 


D  D 

_xau^Q     9i+_ii92L 

R   ax        X  ax       R   ae 


M, 


R  \ap      R      R  a^y      "^  ax 


D 


M  ^     =       xi    a^  +^      +  D      ^ 

x0  R    \de      dx  x6  ax 


D 


M 


ex 


M  /M   _  1  9ii\    +  n       9ZL 

R    lae      R  ae  i        xe  ax 


«x       =     Ox*  -If 


G, 


aw 


%        =     ^    Rr,   .^   - 


R 


86 


(2-11) 
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where 

D 
> 

D 


xe 


E 


Ex     I   , 

''s 

=      E, 

^^3' 

D 

p 

— 

^^oK 

^^60   h^ 

E 

X 

=      E^^h, 

% 

= 

^Oo^ 

E.^h, 

Ox« 

=     Gxe^h, 

G 

X 

= 

k/  GxZq  h 

^/  ""ezo 

h 

1  + 

h^    1 
12R^_ 

» 

I2 

= 

h" 
12 

(2-12) 


The  coefficients  K    and  K„  are  constants  for  adjustment  and  can  be  taken  as 

X  e 


7r/vT2  as  discussed  by  Mirsky  and  Herrmann  (11). 

Substituting  Equations  (2-11)  into  the  Equations  of  Motion,  (2-2)  through  (2-6), 
yields 


2         n       p         F    /  p 


8x 


^    <.x= 


R  \8x 


9x90 


R 


'xe  dxde 


^xe  d^ii) 

2  2 

R     de 


x9  ^    xe\  gfu 

,    ^  R V  9/ 


+   T  ^-!i    +  —   ^-^        —   — 

9x"      r"  a/  "  ^  9^ 


=     h 


Po  "    "^    R    '^ 


R 


(2-13) 


R  \   R 


D, 


R^ 


8w    ^  9^v 

.9^      ^0^ 


D  ^2  E  ^2  2 

e    9  T?    ^  _i^    9  u     ^  9_v 

R=    9/         ^    9^^^  ""^  9x^ 


_^  _xe_   9_2.  +  _x£    9  "     +  _i 
R      ^  2  R     9x961  R 

dX 


r]   +  =- 


1   9w       v_\  ^  „  gfv 
R    90     -r)  3^2 


N 
2 

R     L 


9w    ^  gfv    _^  ,     / 1    9w     ^ 
90         g^a  o\r    90         '^^ 


=     p^  h  V  +  ^  r, 


(2-14) 
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^  '  \       do  J      ^3  do       R   ax        xyax      ^^2 


.R^      R-* 


R 


N 

r" 


ot?  ^  1  /  a  w      av 


^  ^  afw  _  _N,  L+2^  -  ^ 
ax^       r"   V       ^^      a/, 


-f    l^^^ 


l\   +  h     f  1    +  !!^W    +   R    9H    +   Rh     ^1 

ft  y       oy        Rya0  ax      "  o  ax  J 


+  q'(x,  0,  t)     =     ^w  +  p     hw 


(2-15) 


^   ax^        ^  ax^       ^   9^9^       R^  W  ~  ^  do^ 


D  2 

+    xe    a  T? 


R     axa( 


°x^^lf 


Nh 
o 

R 


aw 

ax 


=     Iip  +  ^  u 


(2-16) 


^  (  9liL     1  9w  _  j_  afv\  ^  ^  _af£  ^  ^\o_l_£±      afv ' 
R^ae^      Rae    -R  ^^^^       r   axae        r  ^axao      ^^2 


+  D      ^ 

X0    ^    2 

ax 


0  [^     ,  aw 

■r  r^  ^  ^  -  ^ 


Nh 


R 


2 


^-^^^  If^^--)^ 


It,   +^v 


(2-17) 


Collecting  coefficients  on  deflections  gives 


E     ^^  +  ^  ( G       +  ^  )   -9l   +  T  -^  +  ^    -^ 

^  ax"      r"  \  ^^      rV  a/         ax"      r"   ae" 


9^  h    9M 

^  -  p    h  —    u 

ae  at  J 


(E     +  G     )  -^ — 

^    V         xo'  dxdO 


V   + 


ff<^.-N)^] 


w 


\o    a^ 


2  3  P 

ax        R      do 


1  £. 
«  at^ 


4)    =    0 


(2-18) 
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1  a 


u  + 


—  lE^  +  —   +   N  )    -^   +   (G   ^  +  T)  -^ 


R' 


R 


do 


xe 


R' 


z(^^^M^^^i^(-^)u- 


-    p      h   -^       V 


D^  ^2         D   ^  2  G,        h     N 

e  a     ^    xe  a     +  _£  +  _o 

^3           2                R  2  R                _2 

L    R  ae  ax  r 


I     8=1 


ax 


(2-19) 


^(^.-^>ir 


u   + 


M^. -'i^) -i  ■(•■'*)  t 
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4    k  +  ^   +  N  ^)    -   (G^  +  T)  ^    -  4  (G^  N)  ^ 

R^  \  ^     R^  R  /        X         ax^      R^     ^         ae^ 


^^  ^  ^P 


o'iVi- 


X       R   i  ax 


r4^ 


+  RG,  +  Nh      1  +    „ 
Q  o\  R 


o\l    al 

rjj  a^  h 


q'(x,  0,  t) 


(2-20) 


Kjl    ^i    lill    A 

L  ^   ax^  "   R=   ae^  "  ^  at^J  ""     L 


Nh 


\    X         R 


o  \  _a_ 

ax 


w 


J.  (D     +  D         -^ — 

R  ^  V      xe)  dxde 


TJ        =        0 


(2-21) 


13 


D  2 

xO     d 


n 


Dx 


D.     C.2  ,    /  Nh 

_£  1_   +  1     G     +  -=^ 


1  iC 


_i 

R 


-^  +  G     + 

r"         ' 


Nh 
o 

R 


1   + 


R 


_a_ 


w 


^  u^    xe'  dxde 


] 


2  2  0 

R        80  " 


Nh 
o 

R 


,-Tf  ^^ 


-     I 


at 


0     (2-22) 


Equations  (2-16),  (2-17),  and  (2-18)  express  the  principle  of  linear  momen- 
tum in  the  x,  O,   and  z  directions,  respectively,  while  Equations  (2-21)  and 
(2-22)  express  the  principle  of  angular  momentum  about  an  axis  through  the 

middle  surface  in  the  direction  of  K    and  K  ,  respectively. 

a  X 


These  equations  can  be  written  in  the  form 


C.  u. 


^i' 


ij 


C. 


In  matrix  notation  the  set  has  the  form 


(2-23) 


^11        *^12        ^13        ^14 


^12        ^22        ^23 


25 


^13       ^23       ^33       ^34        ^35 


;         0      c       c       c 

14  34  44  45 


0      c       c       c       c 

25  35  45  55 


u 

0 

V 

0 

w 

= 

q'(x,  e. 

t) 

4' 

0 

V 

0 

(2-24) 
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where  the  coefficients  are  defined  as  follows : 

C  =     (E     +  T)  -^  +  —  (  G       +  -^  +  N  I  -^ 


°     at 


12 


1  a 

R  ^  [^       xe'  dxde 


C  =       :^(E       -   N)    7^ 

13  R  ^     V  '    dx 


C  =     -rf 


D  2  D  2  2 

X   a  xe   a         i    a 


14 


^   ax=       r"  do^  "  ^  at^ 


22 


""^         ax^      R^  V  ^     R^  de~      R~       "     af 


2      -    —     -    Po^    -.2 


c        =    -^ 


23 


25 


R 
D 


■=(.  "  ^  *  °e  *  ^  1,1  *  f 


_a_ 
ae 


^    ax^      R=  ae^       ^       R^       ^  at^ 


33 


-(G^  +  T)  ^   -  -^  (G    +  N)  ^   +  ^      E     +  ^   +  N  ^ 

"^         ax^      R^     ^         ae       R^  \  ^     R^  ^  / 


J.  >•   9    _i      u  a 


34 


-K^'^f^ 


35 


R 


D  /         h 

^.RG^.Nh    (l.f 


d_ 

dd 


44 


D     -9l   +^  ^ 


^  ax^ 


R^     80^ 


G     -  N  ^ 
x  R 


at= 


45 


^     (D    +  D     )  — — 

^  p      xe   dxde 


R 


55  xe 


■-^)-S 


(2-25) 


15 

Axisymmetric  Loading 

Assuming  the  loading  on  the  cylinder  is  axisymmetric,  the  set  of  Equations 
(2-24)  reduces  to  the  following. 

(E     +  T)  ^    -  p     h  ^   +  i-  (E     -  N)  ^   +  ^  ^   -  1   9!i      =     0 
ox  at  9x  9t 

-L/F       \r\  ^       /no-  T-\  9  w   ,   .  aw    ,        ,    a  w 

R  <^.  -  ^)  ^  -  (^x  ""  T)  77  ^  ^  aT  ^  Po  h  -J 

ox  at 

^afuiafu/^   ^N^\|w^D    9!i.i9!i 
^   ax^      ^  at=       \  ^         ^  /  ^^        ^  ax^        at^ 

-  (g^  +   N  ^J^     -      0  (2-26) 


Pressure  Loading  Form 

A  step  input  in  external  pressure  which  is  finite  in  both  magnitude  and  time 
and  which  travels  down  the  length  of  the  cylinder  at  constant  velocity  V  can  be 
represented  in  the  form 

q'(x,  t)     =      -  q{H[Vt  -  x]    -  H[V(t  -  t^)  -  x]  }  (2-27) 

where  H(y)  is  the  Heavi side  step  function,  defined  as 

0,  y  <  0 
H(y)     -     {  (2-28) 

1,  y  >  0 
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This  pressure  loading  is  represented  schematically  in  Figure  2.2 


TTTiTtTTTtTTTTTYT 


-►      V 


R 


-. Vt 


Figure  2.2.    Pressure  Loading 

Nondimensional  Equations 

The  steady  state  solution  will  be  investigated.    Making  the  transformation 
a     =     X  -  Vt  (2-29) 

the  partial  derivatives  may  be  written  in  terms  of  a. 


d 

9x 

d     da 
da   dx 

d 

aa 

d 
8t 

d     da 
da   8t 

=      -V    ^ 

da 

d^ 

=      v^    9^ 

9t  da 

Using  these  relationships  in  Equations  (2-26),  letting 


"  ^  I 


(2-30) 


W    =     ^ 
R 


<t>     = 


R 


(2-31) 
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multiplying  by  R/E    in  the  first  two  equations  and  1/E    in  the  third  gives 


1  + 


jr  \  d"u 

X  /  d0 


2tt  P       h  V  j2 


din   .p 

2  I    E 


d<^' 


N    \  dW 
E  Jd0 


D  ,2, 


R^E      d0^ 
X       ^ 


E     R^    d0^ 

X 


Ejd<^ 


2x    ^  X  1  dfw 

E  E     /     ,  ,2 

X  X  /    dcp 


tVR    dW   ^  PQ  ^  dfw 

E        dd)  E  ,,2 

X        ^  X         d0 


/E,  D, 

\     X        E     R 


^        ^x    ^  / 


2x  ,  ji  ^\  d£ 

E  E       R /    d0 


^  [H(-R0)   -  H(-R<^-  Vt^)] 


D 


d^U 


V^I     d^U 


R^E      d(f 
x 


R^E 


d(^' 


fx    +  Ji  l^o  \  dW 
E  E       R    I    d(^ 

X  X  / 


D 


R^E 


df£ 
d(^^ 


E    R^   d0^ 
X  ^ 


X    ^     N      o  \    ,     _     rt 

r-   ^  E"    R-j'^     -     ° 
X  X         / 


(2-32) 
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Now  the  following  dimensionless  ratios  are  defined. 
T 


X 

E 

X 

F 

N 
E 

X 

= 

P 

E 

X 

= 

% 

E 

u 

E 

X 

= 

\' 

G 

X 

E 
E 

h 
o 

R 

= 

r 
o 

x" 

= 

'^o 

R 

E 

X 

r 

= 

P    h 
^o 

^0 

h 
R 

I 

o 

= 

I 

12R=  P„ 

h= 

^o«^ 

12R^ 

2 

= 

X    \         E     R^  /         R^  E  ^ 


(2-33) 

The  dimensionless  velocity  parameter  is  denoted  by  X  and  r  refers  to 
the  inertia.    A  distinction  will  be  made  in  the  axial  and  radial  inertia  terms  for 
later  investigation,  r    denoting  axial  and  r  denoting  the  radial  inertia  terms. 
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Using  the  dimensionless  ratios  given  in  Equations  (2-33),  the  equations  of  motion 
now  take  the  form 

(1  +  F  -  r    x^)  ^   +  (E     -  p)  1^  +  (D     -  I    A^)  ^     =     0  (2-34a) 


(E     _p)^   _(G+F-rA")^^   -.A^   ME     4-Pr)W 

-  (G  +  P  r^)  ^   =    -q^[H(-R0)  -  H(-R(^  -  Vt^)]  (2-34b) 


(D     -  I    x'")  ^   -  (G  +  P  r  )  ^  MD    -  I    A^)  ^ 
^        °  d0^  o'    d0         ^    1        o       '  ^^2 

-  (G  +  P  r  )i^     =     0  (2-34C) 


CHAPTER  ni 
TRANSFORMATION  OF  EQUATIONS 


Equations  (2-34)  will  be  transformed  using  the  Fourier  transform 

00 

F[f(0)]      -     f  (s)     =      J   f(cp)  e"'^'^  d0  (3-1) 


The  inverse  transform  is 


00 

■"■■■  ds  (3-2) 


F  V'(s)]     =     i{<P)     =     27    /  ^"^^)  ^'^"^ 


where  i  =  yf-T . 

Assuming  all  the  derivatives  of  f(<^)  through  order  (r-1)  vanish  sls  (p  -*  ±«> 
the  transforms  of  derivatives  of  f((/))  are  given  by 

{\s)     =    (is)^f(s)  (3-3) 


so 


I    dcf.^  J 


is  f (s)  (3-4) 


F  ^^        -     -s^f(s)  (3-5) 
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Applying  this  to  Equations  (2-34)  gives 

(1  +  F  -  r  \^)  s^  U  (s)   -  i(E     -  P)  s  W  (s)  +  (D     -  I    X^)  s^  4(s) 


i(E    -  P)  s  U  (s)  +  (G  +  F  -  r  X^)  s^  W(s)  -  i  e  X  s  W(s) 


+  (E^  +  P  r^)  W  (s)  -  i(G  +  P  r^)  s  ip  (s) 


^o  ,     isd      ,,  ^ 
IS  '  ' 


(D^  -  Iq  ^^)  s^  U(s)  +  i(G  +  P  r^)  s  W  (s)  +  (D^  -  I^  X^)  s^  4(s) 


+  (G  +  P  r^)  i^  (s)     =     0 


(3-6) 


Collecting  coefficients  on  like  displacements  allows  this  set  of  equations  to 
be  written  in  the  following  matrix  notation. 


[(l  +  F-r^X^)s^]      [-i(E^-P)s]        [(D^-I^X^)s^] 


[(G  +  F-rX^)s^ 


[i(E^-P)s]         -i€Xs 


+  (E    +Pr    )] 
^    o  o 


2,    2 


[-i(G  +  Pr^)s) 


[(D^-I^X^)s^ 


[(D^-I^X)sl       [i(G.Pr^)s]         ,(G.Pr)] 


U(s) 


W(s) 


i^(s) 


0 


o  ,    isd    , , 
•:—     e       -1 

1  s  ^  ' 


(3-7) 


*See  Appendix  A  for  the  derivation  of  the  transform  of  the  forcing  function. 
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This  set  of  equations  can  now  be  solved  for  U,  W,  and  ijj.    This  work  is  carried 
out  in  Appendix  B  and  the  results  are 

^^-^     =     ^^  ~  ^ — ^  [(D     -  I    A^)  (G  +  P  r     +  E     -  P)s^  +  (G  +  P  r  )  (E    -  P)] 
q  2  TT,   ,     ^    1       o  ^  o         1  ^  o'  '    1        ' 

^o  s     D(s) 


^^^    =     ^^  "     — ^  [(D     -  I    X^)c  s^  +  (G  +  P  r   )  (1  +  F  -  r    A^)] 
q  .     — ,   ,      ^^    1       o       '  ^  o' ^  1      '^ 

^o  isD(s) 


^^      =     -  ^ — — ^  [(E     -  P)(D    -  I    A^)   +  (G  +  P  r   )(1  +  F  -  r    A^)] 

^o  D (s) 

(3-8) 

where 

—  4  3  2 

D(s)     =     c     s     +ic     s     +c     s     +ic     s  +  c  (3-9) 

and 

c       =     (G  +  P  r  )  [(1  +  F  -  r    A^)  (E     +  P  r   )   -  (E     -  P)^] 

0  o  1^0  01 

c       =     -  e  A(G  +  P  r   )(1  +  F  -  r    A^) 

1  ^  o^  1       ' 

c       =     (D    -  1    A^)  {(E     +  P  r  )c   -  (E     -  P)  [(E     -  P)  +  2(G  +  P  r  )]  } 

2  1         o  "^      O  O  1  1  o  ^ 

+  (G  +  P  r  )  (1  +  F  -  r,  A^)  (F  -  P  r     -  r  A^) 

c       =     -  e  A(D    -  I    A^)c 

3  ^10 

c       =     (G  +  F  -  r  A^)(D^  -  I    A^)c 

4  ^  '^10 


c       =     1  +  F-D^-r^A^  +  I^A^  ^^_^^^ 

The  displacements  are  found  by  inverting  Equations  (3-8)  using  transforma- 
tion (3-2).    In  order  to  evaluate  these  integrals  the  roots  of  the  characteristic 
equation  D(s)  =  0  must  be  determined. 


CHAPTER  rV 
INVESTIGATION  OF  THE  CRITICAL  VELOCITIES 

It  is  known  that  the  inverse  of  the  deflections  [given  by  integral  (3-2) ,  where 
f  (s)  represents  the  deflection  expressions  (3-8)]  does  not  exist  when  there  are 
repeated  roots  of  the  characteristic  equation 

D(s)     =     0 
on  the  real  axis.    This  can  occur  when  there  is  no  external  damping  and  corres- 
ponds to  a  resonant  condition  as  discussed  by  Jones  and  Bhuta  (2).    There  are 
specific  load  velocities  corresponding  to  these  points  and  they  will  be  referred  to  as 
critical  velocities.     The  condition  which  must  be  satisfied  in  order  to  have  re- 
peated roots  is  that  the  discriminant  of  the  undamped  characteristic  equation 

cs'^+cs^+c=0  (4-1) 

4  2  0 

must  be  zero.    The  discriminant  of  this  equation  is  determined  in  Appendix  C  as 
A     =     16c    c    (c^    -    4c    c  )^  (4-2) 

O      4  ^     2  O      4' 

Therefore,  the  three  conditions  which  will  make  the  discriminant  zero  are 

c^    -    4c    c       =0  (4-3) 

2  O     4 

c^     =     0  (4-4) 

4 

C        =     0  (4-5) 

o 

Substituting  the  required  coefficients  from  (3-10)  into  (4-3)  gives* 

C   X^    +   C\^    +   C   X^    +   CX^    +   C       =0  (4-6) 

8  6  4  2  O 


*See  Appendix  D  for  the  detailed  calculations. 
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24 

where 

C       =     e^    -    e 

8  3  14 


C       =     2e    e      -    e 

6  8     9  13 


c 

4 

= 

e^    +   2e  e      - 

a             7  9 

\^ 

c 

2 

= 

^^7%     -     ^1 

c 

o 

= 

e^    -    e 

7                 10 

The  e.  are  defined  as 


8         =De     +eff 

7  15  3    O    2 


6         =Dee      -le      -e„f^ 

8  104  05  35 


e         =err     -lee 

9  3         1  O     O     4 


e         =     4e    e    f   D   f 

lO  3     6    3      11 


(4-7) 


e         =     4e    [e,(ff    -rDf)    -erfDf] 

11  3^6^39  11'  01311^ 

e         =     -4e    [e,  (rf    +  f  I   e  )    +   e   r,  (f  f    -  rD,  f )] 

12  3^6^        9  30     4'  0139  ll'^ 

e         =    4e    [e    rl  e     +   e   r    (rf     +f  I  e  )] 

13  3         6        O     4  O     1  ^        9  3    0     4'-' 

e,  ^     =     -4e^e^r    rl   e^  (4-8) 

14  30104  ^  ' 

and  the  f.  coefficients  are  defined  in  Appendix  D. 

If  the  axial  and  rotatory  inertia  are  neglected  the  equation  for  determining 
the  location  of  the  repeated  roots  reduces  to 

CA^+CA^+C      =0  (4-9) 

4  2  O 


where 


t;  2  2^ 

C      =     e    r   I 

4  3  O 


c 


2 


=     2re    [2e,  Df     -f(De+eff)] 

3'         611  0^15  302' 


5      =(De+eff)-4efDfe,  (4-10) 

O  ^1530    2'  33116  ^  ' 
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The  solution  of  Equation  (4-9)  gives  two  roots  (X'  )^  and  (X'  )^  which  are 

1  2 

approximations  for  X^  and  X   .    These  are  the  square  of  the  two  lowest  load 
velocities  at  which  repeated  roots  occur. 

The  second  condition  which  gives  repeated  roots  is  that  given  by  Equation 
(4-4).     Substituting  for  the  c^  coefficient  gives 

(G+F-rX^)  (D^  -I^X^)  [l+F-D^  +  (I^-r^)X^]      =     0  (4-11) 

This  can  be  written 

(f^  -  rX^)  (D^  -  I^X^)  (f^  +  e^X^)     =     0  (4-12) 

2 

Expanding  this  and  collecting  coefficients  on  X    yields 

[f3D^-(rD^-.f3yx^    +    I^rX^]   (f^-.e^X^)     =     0 

+    [f   rl      -    e    (rD    +f   I  )]X^    +   e    I   rX^        =     0 

lO  4^1  3    0"  40 


Finally 


C'  X^   +   C'  X*    +    C'  X^    +    C'      =0  (4-13) 

6  4  2  O  ^  ' 


w^here 


C'  =     e    I   r 

6  4    O 

C'  =     f   rl     -    e    (rD    +f   IJ 

4  10  4^130' 

C;  =     f^f3D^  (4-14) 
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If  rotatory  and  axial  inertia  are  neglected  Equation  (4-13)  reduces  to  give 
one  root 

C'  f  p 

The  third  condition  (c    =  0)  from  Equation  (4-5)  is  now  investigated.    Substituting 
for  c    gives 

(1  +  F-  r^x"^)  (E^  +  Pr^)    -    (E^  -  P)^     =     0  (4-16) 


or 


(f    -  r   A.^)  e       =     e  (4-17) 

^  o        1      '    o  1  ^         ' 


thus 


A 

e    f    -  e 

O    O           1 

e    r 

O     1 

K 

- 

^6 

e    r 

O     1 

(4-18) 

Therefore,  Equations  (4-6),  (4-13)  and  (4-18)  can  be  solved  to  give  eight 
roots  in  }?  which  satisfy  the  conditions  for  repeated  roots.    Equation  (4-6)  gives 
four  values  of  \^  and  these  values  are  labeled  A^,  a|,  X^,  and  a|.     The  condi- 
tion which  lead  to  these  roots  was  that 

c^    -    4c    c       =0 

2  O      4 

Solving  Equation  (4-1)  directly  gives 


-  c    ±  ,/c^  -  4c    c 
s^     =  ^      ^    ^ 2_^  (4-19) 


If  the  radical  in  (4-19)  is  zero,  the  roots  are 


^.= '  y-  ^  •       ^.. '  -  V  -  ^       <^-^°' 
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The  repeated  roots  are  either  real  or  imaginary  depending  upon  the  sign  of 
the  coefficients. 

Equation  (4-13)  gives  three  more  roots  which  will  be  labeled  A. 3,  X^,  and 

2 
A.^.    These  roots  come  from  the  statement  that  c^  =  0  and  if  c    /  0  the  repeated 

roots  in  this  case  will  occur  at  an  infinite  value. 

Equation  (4-18)  adds  one  more  root,  making  a  total  of  eight.    This  root  is 

2 
designated  Xg.    The  condition  leading  to  this  root  was  that  c    =0.     From  Equa- 
tion (4-1)  it  is  observable  that  the  characteristic  equation  becomes 

s^ic^s^   +   c^)     =0  (4-21) 

which  shows  a  repeated  root  at  the  origin. 

A  computer  program  was  written  for  the  solution  of  these  equations.    A 
simplified  flow  diagram  of  the  program  is  shown  in  Figure  4. 1  and  the  details 
of  the  program,  named  VCRIT,  are  presented  in  Appendix  E. 

The  results  of  a  parametric  study  using  the  computer  program  VCRIT  are 
presented  in  Figures  4.  2  through  4.  8.    The  locus  of  each  of  the  roots  A  ^(1  =  1,   ...8) 
is  shown  on  these  figures.    These  curves  are  the  boundaries  which  separate 
these  plots  into  distinct  regions  which  are  labeled  as  Regions  I  through  Vin.    In 
each  of  the  regions  the  roots  of  the  undamped  characteristic  equation  have  a  par- 
ticular form  as  noted  on  the  figures  and  as  listed  in  Table  4-1. 

Table  4-1 
Correlation  of  Root  Type  with  Region  Numbers  for  Figures  4.  2  through  4.  8 


Region 

Form  of  Root 

I.  m,  VI 

±ibi,  ±ib2 

II 

a  ±  ib,  -  a  ±  ib 

IV,  VII 

±ai,  ia^ 

V,  vni 

±a,  ±ib 

a,  1 

0  are  real 
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Insert  Data  in  Separate  File 

z        E,, 


h     ^Oq     Gxz„      E^,^ 


R  '   Ev     '    Ey      '   Ey 

Aq  Xq  Aq 


P,   F 


Start  VCRIT 


Yes 


Read  and  Print  Data 


Calculate 
Coefficients 


Solve  Fourth  Order 
Equation  in  A^ 


Print  Roots 


Solve  Third  Order 
Equation  in  X^ 


Print  Roots 


Solve  First  Order 
Equation  for  Last  Root 


Set  Axial  and 

Rotatory  Inertia 

to  Zero 


Repeat  Solutions  Using 
Simplified  Equations 


More  Data? 


No 


Stop 


Figure  4. 1.    Flow  Diagram  of  Computer  Program  VCRIT  which  Determines 
Load  Velocities  at  which  Repeated  Roots  Occur 
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10.000 


1.000 


100 


0.1 


0.01    ^ 


0.001  ^"'•■•■'•''•' 
0.00001 


Pure 

Imaginary 

Roots 


1.0.  t     =    0 


TT-ii  =    0.3.  -^    =    0.35 


P    =    0.004.   F  =  0.002 


2  Real, 

2  Imaginary 

Roots 


0.0001 


0.001 


0.01 


0.1 


h/R 


Figure  4.2.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  the  Thickness  to  Radius  Ratio  Including  Prestress 
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10,000 


1000 


100 


0.1 


=    0.35 


0.00001  .0001  .001  .01  0.1 

h/R 


Pure 

Imaginary 

Roots 

2  Real, 

2  Imaginary 

Roots 


Figure  4.3.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  the  Thickness  to  Radius  Ratio  with  No  Prestress 
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1000 


Pure 

Imaginary 

Roots 


2  Real, 

2  Imaginary 

Roots 


o.iot^ 


0.40 


0.80 


1 .: 


l.G 


2.0 


Ee^Ex^ 


Figure  4.4.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  Eg   /Ex„ 
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Real  Roots 
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Imaginary. 
Roots 


.00001 
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Figure  4.5.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  G^zq/^Xq 
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VII Real  Roots 
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0.0001 


0.0001 


0.001 


0.01 
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1.0 


10 


Figure  4.6.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  the  E^;  /Ejj^ 
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0.1 


liegion  I 


j       Pure 
"  Imaginary 
Roots 
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'■^■^  ■-^^^^^-' 


-0.002 


0  0.002 
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0.004 


Figure  4,7.     Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  the  Circumferential  Prestress 
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1.0 


0.1 


-0.002  -0.001 
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Figure  4.8.    Classification  of  Roots  of  the  Undamped  Characteristic  Equation  for 
Variations  in  the  Axial  Prestress 
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The  lower  three  values  of  X.  are  approximated  by  Expressions  (4-9)  and 
(4-15)  by  neglecting  the  axial  and  rotatory  inertia.    The  approximations  are 
quite  good  over  the  range  of  parameters  studied. 

The  critical  load  speeds  are  denoted  as  those  speeds  at  which  the  displace- 
ments become  unbounded  for  an  undamped  system.    This  corresponds  to  the 
load  speeds  which  produce  a  double  root  on  the  real  axis  as  will  be  shown  later. 
It  is  instructive  to  follow  the  path  of  the  roots  of  the  undamped  characteristic 
equation  in  the  complex  plane  as  the  load  speed  increases.    As  an  example, 
the  roots  will  be  traced  for  material  properties  corresponding  to  an  isotropic 
shell  with  h/R  =  0.  001  and  positive  axial  and  circumferential  prestress. 

Following  the  vertical  line  for  h/R  =  0.001  in  Figure  4.2  for  increasing  X^ 
will  give  a  path  crossing  all  of  the  boundaries  separating  different  types  of 
roots.    Starting  at  the  low  load  speed,  the  roots  are  all  on  the  imaginary  axis. 
(This  would  not  be  the  case  if  prestress  were  not  included,  as  shown  by  Fig- 
ure 4.3.)    As  a  means  of  tracing  the  location  of  the  roots.  Figure  4.9  is  util- 
ized which  shows  the  complex  plane.    The  roots  in  Region  I  appear  on  the  plane 
on  the  imaginary  axis  and  these  particular  roots  are  designated  s       ,  s       , 
s„  ,  ,  and  S^   ,  .    The  nomenclature  s.   .  denotes  the  i^^  root  location  and  j  indi- 

3,1'  4,1  l,j  J 

cates  the  relative  position  of  the  roots.    For  instance  s     ^  denotes  the  position 
of  the  first  root  and  at  that  time  (load  speed)  the  other  roots  are  located  at 
points  s^  ^ ,  s„     ,  and  s,    , .    The  arrows  indicate  the  direction  in  which  the 

■^  2,4'       3,4  4,4 

root  is  moving  for  an  increasing  load  speed .  At  the  first  speed  two  pairs  of  roots  are 
moving  toward  one  another  on  the  positive  and  negative  parts  of  the  imaginary 
axis.    They  meet,  and  the  first  repeated  root  location  is  established  which  cor- 
responds to  X    in  Figure  4.2.    As  the  load  speed  increases  Region  II  is  entered. 
The  roots  are  complex  as  can  be  seen  in  Figure  4.9.    Next,  the  complex  roots 
approach  one  another  in  pairs  on  the  negative  and  positive  real  axis.     This  gives 
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the  first  repeated  roots  on  the  real  axis  and  this  load  speed  is  designated  as  the 
first  critical  load  speed,  ^ipr>. 


Complex  8  Plane 


8.   .  =  Root  1  at  load  speed  j 


Figure  4.9.     Path  of  the  Roots  of  the  Undamped  Characteristic  Equation 
for  Increasing  Load  Speed  (Beginning  in  Region  I) 
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Now  the  roots  separate  and  go  in  opposite  directions  along  the  real  axis ,   the 
larger  roots  eventually  becoming  unbounded.    This  speed  corresponds  to  the 
boundary  line  between  Regions  IV  and  V.     The  condition  causing  this  occurrence 
is  that  c^ -*•  0.    Repeated  roots  occur  at  infinity  and  the  second  critical  load 
speed  has  been  determined.    As  the  load  speed  increases  further,  the  large  roots 
come  in  from  +i<x'  along  the  imaginary  axis.    At  location  5  the  two  large  roots 
are  imaginary  and  the  other  two  are  still  on  the  real  axis  approaching  the  origin. 
This  corresponds  to  Region  V  in  Figure  4.2. 

The  locations  of  the  roots  are.now  transferred  to  Figure  4.10  to  avoid  undue 
complication  of  the  picture.    Two  roots  meet  at  the  origin  while  the  other  two 
are  yet  large  and  imaginary.    This  double  root  corresponds  to  A^  in  Figure  4.2 
and  is  the  third  critical  load  speed.    An  increase  in  load  speed  now  causes  the 
roots  to  become  all  imaginary  which  is  in  Region  VI.     The  last  four  values  of 

2  2 

X.  (1  =  5,6,7,8)  are  very  close  together.    In  fact  X    is  approximately  1000,  and 

2  2  2 

X^,  X^,  and  X^  are  almost  nondistinguishable  at  1002.    Increasing  the  load  speed 
past  X^  causes  the  large  imaginary  roots  to  proceed  to  an  unbounded  imaginary 
value  and  reappear  on  the  real  axis  so  that  the  roots  are  now  two  real  and  two 
imaginary.    Also  X    must  be  a  critical  speed  because  Condition  (4-4)  is  satisfied 
at  load  speed  X        This  region  is  not  distinguishable  on  Figures  4.  2  through  4.  8 
and  is  not  given  a  number. 

Further  increase  in  load  speed  moves  the  imaginary  roots  out  and  eventually 
they  reappear  on  the  real  axis  with  the  other  pair,  there  to  remain.    This  posi- 
tion is  indicated  by  s.         (1  =  1,2,3,4).    It  is  interesting  to  note  that  X^ 

1 , lO  6,7,8 

also  must  correspond  to  a  critical  load  speed,  X^^-^. 

The  lowest  critical  load  speed  for  an  isotropic  shell  is  therefore  X^pp  =  X^. 
Figure  4.5  shows  that  for  a  material  with  a  very  low  shear  modulus,  the  first 
critical  speed  becomes  X^  =  X^^ri  which  corresponds  to  the  shear  wave  speed. 
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Figure  4. 10.    Path  of  the  Roots  of  the  Undamped  Characteristic  Equation  for 
Increasing  Load  Speed  (Beginning  in  Region  V) 
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When  this  load  speed  goes  to  zero,  the  critical  small  deflection  buckling  load 
has  been  reached.     From  Equation  (4-15)  this  load  is  seen  to  be  approximately 

F     =     -G  (4-22) 

This  corresponds  to  the  second  critical  load  for  engineering  materials  with 
a  more  realistic  value  for  shear  modulus. 

If  prestress  is  not  included,  Region  I  disappears  as  shown  by  comparing 
Figures  4.2  and  4.3.    Decreasing  the  tangential  modulus  drastically  or  increas- 
ing the  normal  modulus  results  in  the  same  effect  as  shown  by  Figures  4 . 4  and 
4.  6.    Circumferential  prestress  has  essentially  no  effect  on  the  critical  load 
speeds  but  axial  prestress  has  a  pronounced  effect  as  shown  by  Figures  4 .  7  and  4.8. 

It  is  instructive  to  make  a  comparison  of  these  results  with  those  obtained 
in  Reference  (2).    Threecritical  velocities  were  derived  there,  which  are  given  as 

p^R[3(l-M)]2         6p^R^(l-M) 


VcRs     =     ^  (4-24) 


VCR3     =    f  (4-25) 

'^o 

for  an  isotropic  material.    Since,  in  the  nomenclature  used  in  the  present  work, 

V^pR(l-M^) 
^       =    -^Eh (4-26) 


2 

The  corresponding  expressions  in  terms  of  A.    are 


^CRs     =     f  (4-28) 

^CR.     =     #(1-M^)  (4-29) 
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Now  taking  /i  =  .  3  it  is  found  that 

^CR        =     0.55075    -    0.15^  (4-30) 

^CR^     =     f  (4-31) 

^CR3     =     0.9lf 

2  h 

From  Figure  4.  3  (which  is  a  plot  of  \.  versus  -^  for  an  isotropic  material 

with  ^  =  .  3)  for  the  case  of  zero  axial  and  circumferential  initial  stress,  A(^j^  , 

2  2  2  2 

'^CRc»  ^"^  ^CR-^  ^^^^  Reference  2  agree  extremely  well  with  \   ,  X  ,  and 

A^  respectively.    The  two  velocities,  A    and  A    ,  arise  in  the  present  results  be- 
cause of  initial  prestress  considerations  and  shear  deformation,  respectively, 
which  were  not  included  in  the  referenced  results.    A    corresponds  to  the  dili- 
tational  wave  speed,  A^  corresponds  to  the  shear  wave  speed,  A^  corresponds 
to  the  bar  wave  speed,  and  A     ,       ^  corresponds  to  the  plate  wave  speed. 

5 , o , y ,8 


CHAPTER  V 
SOLUTION  FOR  DISPLACEMENTS 

General  Solution 

Transformed  Displacements 

From  Equations  (3-8)  the  transformed  displacements  can  be  written  as 

2 

TT/  \            1        isd     c'   s    +  c' 
U(s)      ^     1  -  e  _5 6 

%  s^  D(s) 

m7/<,\            1       „isd     c'   s    +  c' 
W(s)     ^     1  -  e  _7 8 


^o  D(s) 


q  (1  -  e^-)  -^  (5-1) 


^o  D(s) 

where 


c'      =     f  (e    +  e  ) 

5  7^    3  l' 


c'  =     e„  e 

6  3      1 

c'  =     f  (f    +  e    x^) 

7  7^    1  4  ' 

C  =     e  (f^  -  r    X^) 

8  3     O  1  ' 

c'  =     e    f    +  e  (f    -  r    A^) 

9  17  3^  O  1         ' 

f  =     D     -  I    ;^^ 
7  1        o  (5-2) 
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Inverse  Transformation  of  the  Rotation 

Using  the  inverse  transformation  given  by  Equation  (3-2)  the  rotation  is 
defined  as 

^     —  I    /  M^  e'*^  d.  ,5-3, 

%  ^"^     ^  D(s) 

By  partial  fraction  expansion  this  can  be  put  in  the  form* 


where 


Q^k  "   T"; ; "   ~4 ^^^ 

n  (\  -  ^j 


|At^(«)i 


s=s,  --km 

k  m=l 


c\ 


(5-4) 


%  =  r  (5-5) 

4 

Defining 

CO!,         /«    /I        isd,      .  , 

—  CO 

then 

4 
^       ^      y;     F,(c^)  (5-7) 


*See  Appendix  F  for  the  derivation  of  a  sample  partial  fraction  expansion 
expression. 
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Letting 


1        isd 

f,  (s       =     e  ^  (5-8) 

k 

the  integral  to  be  evaluated  is 

00 

f  ys)  ds 

_  00 

where  fi^(s)  is  an  analytic  function  except  at  the  simple  pole  s  =  s.  .    Defining  the 
complex  root  in  general  to  be 

\     =     ^k  ^  ^  \  • 
the  Cauchy  integral  theorem  is  used  to  evaluate  this  integral . 


Im(s) 

/v 

^^'^ 

l\ 

T— > 

^k 

'] 

^\^ 

Re( 

Figure  5.1.    Contour  Integration  Path  for  Evaluating  Rotation  Integral 


Assuming  first  that  b,    >  0  the  Cauchy  integral  theorem  gives 


i 


f,  (s)  ds     =     27r  i  2  Residues 


(5-9) 
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The  integral  around  the  closed  path  shown  in  Figure  5.1  is 

R 

<k     f|^(s)  ds     =      y    fj,(s)  ds  +     J     fj^(s)  ds  (5-10) 

From  Reference  (12)  it  is  shown  that  if  f(s)  -*  0  uniformly  as  R  -*  «=  then 

lim        J    i{s)  e^^^  ds     -     0,  (^  >  0) 

R^co     c^ 

Therefore  the  first  integral  on  the  right  in  Equation  (5-10)  goes  to  zero  as 
R  -►  00  and 


(|)      fj^(s)  ds     =        /       f  (s)  ds     =     2ir  i  Res(s.  ) 


—  00 


i0(a^+ib^) 
=     27r  i    e  ,  (j)  >  0 


liril- 


id(a.+ib,)-l    i(^(a  +ib  ) 
+  27ri|-e         "       ''e         ^       ^  ,        0  +  d>O 


Therefore 

-bj^  <^+i  aj^  cp  -h^(d+(t>)+i  a^(d+(P) 


r  -b.  <?)+ia.  0  -b  (d+(/))+ia.(d+(^)  -| 

F^     =     -iCg  cv^    e     *"  •"      H(0)-e     "^  "^  H(d  +  0) 


H(d  +  0)       (5-11) 


where 

\   >    0 

H     =    Heaviside  step  function. 
Similarly,  when  b,   <  0  the  integration  path  is  in  the  lower  half  plane  and  the 
result  is 

j   -h.(t>+ia,(t>  -b  (d+0)+ia,(d+0) 

^k    "     ^''s^kr  H(-<^)  -  e     ""  ""  H[-(d+0)l(, 

bj^  <  0         (5-12) 
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The  general  expression  for  F,   can  be  written  in  the  form 


Fj^  =    -  i  sgn(b^)  c^a^\e 


-h^(d+ct>)+isi^(d+(P) 


H[sgn(bj^)0] 


H[sgn(bj^)(d+(^)] 


where  H  is  the  Heaviside  step  function  and 

1. 


sgn(bj^)     = 


b,    >  0 
k 


-1,         b,    <  0 
k 


(5-13) 


(5-14) 


The  rotation  is  given  by  Equation  (5-7). 
4 


k-l 


Inverse  Transformation  of  the  Radial  Deflection 

The  radial  deflection  is  obtained  by  inverting  W  to  give 


q^ 


—    f  - 

Ini    J 


00  2 

isd  c'   s    +  c' 

-  e  7              e      i<ps    , 

e  ^     ds 


D(s) 


(5-15) 


The  term  (c'   s    +  c' )/D(s)  can  be  expanded  in  partial  fractions  to  give 


WM    = 


iTT  i      J 


-  e 


isd 


4         ^, 


s  -  s, 


i0s    , 
e  ^     ds 


k=l 


where 


(5-16) 


c     s,      +  c        * 

7     k  e 


'k  4 

n  K-^j 

m=l 


k^  m 


(5-17) 


*See  Appendix  E  for  this  expansion. 
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Defining 


w 


k^^'  27r  i      7       s(s  -  s,  )     ^         d 


(5-18) 


The  deflection  is  then  given  by  the  sum 
k=l 


(5-19) 


By  contour  integration  in  the  upper  and  lower  half  planes ,  the  integral 
Equation  (5-18)  can  be  evaluated.    Breaking  the  integral  into  two  parts  gi 


in 
ves 


F,  (0) 


k 

2-n  i 


I 


,i0s 


s(s  -  Sj^) 


ds 


/ 


gi(0+d)s 

s(s  -  Sj^) 


ds 


(5-20) 


The  contour  shown  in  Figure  5.2  is  used  to  evaluate  these  integrals  for  0  >  0, 


b,    >  0, 

k 


Im(s) 

/ 

/^ 

^"Ncr 

t 

XS" 

^  P         k     ) 

-R 

-* 

P 

<-               R 

Figure  5.2.    Integration  Contour  for  Radial  Deflection  for  d)  >  0    b    >  0 

'    k 

From  Cauchy's  integral  theorem 
/cw^^)^^     -     /wy«)ds+    f     ^f^(s)ds^    /wV«)^« 

Cr  -«  c^ 


+  y    ^ys)ds     =     27ril  Res 
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From  Jordan's  Lemma  (12)  it  can  be  seen  that  the  integral  on  contour 


Co  ^OasR-*«>,  i.e., 


lim 


/ 


i(ps 


R->°o    C 


s{s  -  s^) 


ds 


0,  0  >  0 


R 


and 


lim 


/ 


s(s  -  s.  ) 


ds     ->     0,  0+d  >  0 


R->oo    c^ 


Also 


lim 

£-0 


ZP 


R 


P       /        f,  (s)  ds 


which  is  the  Cauchy  principal  value  of  the  improper  integral .    If  the  integral 
exists,  then  this  is  the  correct  value  for  the  integral  and  the  symbol  P  can  be 
dropped.    Therefore 


/ 


i0s 


ds  +  lim 


/ 


i0s 


P 


ds     =     27r  i  Res(s,  ),       (p  >  0    (5-21) 


The  second  integral  in  Equation  (5-21)  can  be  evaluated  as 


lim  / 

p->o     ;; 


i0s 
s(s-sj^)    d^     =     cviRes(O) 


Since  the  integration  is  clockwise 

a     =     -IT 
and 


Res(O)     = 


icps 


s  -  s 


k  /  s=0 
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so 


/i0s 
-^ 
s(s  -  s, 


7r  1 


lim  /      -; r-     ds       = 

r,      J     s(s  -  s,  )  s, 

p-*0     ^        ^  k'  k 


Integ 

ral 

(5-21) 

now 

may 

be  written 

00 

/ 

—  00 

gi^^s 

ds 

=     _ 

7ri 

^  27r  i 

e 

s(s 

-\) 

\ 

where 

0    > 

0 

\    > 

0 

(5-22) 


Similarly 


TTi     ,    27rie^^^^' 

(^  +  d    >    0 

\    >    0 

(5-23) 


Now  the  case  is  investigated  where  b,    >  0 ,  0  <  0 .    The  contour  shown  below 
is  used  for  this  case. 


Re(s) 


Figure  5.3.    Integration  Contour  for  Evaluating  Radial  Deflection  for  0  <  0, 

b,    >  0 
k 
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The  integrals  have  the  values 


/ 


i(ps 


ds  +  lim        / 


i(ps 


_e 

s(s  -  s.  )  ""    ■    ""'       J     s(s  -  s,  ) 


ds     =     0,  (p  <  0 


or 


/ 


^^f^ds     =     -.iRes(O)     =     ^,  c/xO 

b,    >  0 


(5-24) 


Similarly 


/ 


b,>0 


(5-25) 


Equations  (5-22),  (5-23),  (5-24),  and  (5-25)  can  be  combined  to  give  the  solution 


for      FyXcp)  when  b,    >  0  as 


wV*)  =  Irr  1  (- 1;;  ^  ^  »'*'']"<*)  ^  f;;  «<-*> 


;  T  *  ^  ;<*^>M  H(^,d)  + 


|1  H(-*-d)] 


b^>0 


or  finally 


F,  (0)     =     — 
w    k^^'  s, 


/    1  i<^s    \ 

-  I   +  e       '^     H{<P)  +  ^  H(-</)) 


1         H4>+<i)%\  1 

^  +  e  '^     H(04ti)   +1  H(-(|)-d) 


(5-26) 


where  b,    >  0. 
k 
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The  case  where  b,    <  0  is  now  investigated.    Integrating  in  the  lower  half 
plane  it  is  found  that 

i<^s. 


/ 


id)S  /•  id)s 

ds  +     /     TTZ ^T  ds     =     2-ni  — 


s(s  -  s,  )  J     s(s  -  s,  ) 


or 


/ 


"»         .  .  i0s, 

10s  ^  k 

^  ds    =  -27r  i  ^- +  T^  ,  0  <  0 


s(s  -  Sj^) 


s,  s,    ' 

k  k 


bj^<0 


It  is  evident,  therefore,  that  the  general  expression  for       ^\S4>)  can  be  written 
in  the  form 


sgn(bj^)  /3j^ 


w    k^^'  s 


1  ,     ''^\ 
2^' 


H[-sgn(b  )0] 
H[sgn(b^)0]    +  ^ 


r  i(0^)s  -1 

l^-i  +  e  ''J  H[sgn(bj^)(0+d)] 


+  -  H[-sgn(b^)((^4d)] 


Inverse  Transformation  of  the  Axial  Deflection 


The  axial  deflection  in  the  s-plane  can  be  written  as 


(5-27) 


^      s^D(s) 


(5-28) 


^o  =        D(s) 

From  Equation  (5-1)  it  can  be  seen  that  the  first  term  in  Equation  (5-28)  can  be 
written  in  terms  of  ip .    Thus 


q^  c'      q^  2Tr      J        ^2  j^ 


isd.      . , 
e       )     10s   J 
'  e  ^     ds 


(5-29) 


9         'O 


s     D(s) 
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The  integral  in  Equation  (5-29)  can  be  written  as 

4 


2tt      J 


isd,     .     , 


k=l 


a, 


s  -  s. 


ds 


where 

c 

^6      =      C 

Now  defining 


F,  (0) 


c     a 

o 


l^k    r  ii_,^i!ii^,. 


S    (S   -  Sj^) 


then 


(5-30) 


]M1 


_!^  M^  +y-     F,(0) 


q 
9       ^o 


k=l 


The  integrand  in  Equation  (5-30)  can  be  written  in  two  parts  as 


F,  (<^) 


6     k 

27r 


/ 


i0s 


ds 


s    (s  -  Sj^) 


/ 


i(0+d)s 


■ds 


(5-31) 


(5-32) 


s    (s  -  Sj^) 

The  contours  shown  in  Figures  5.2  and  5.3  are  used  to  evaluate  the  integrals 
in  Equation  (5-32)  for  the  case  b,    >  0.    Define 


i0s 


i{(p+d)s 


f,  (s) 
u  k^  ' 


s    (s  -  Sj^) 


f;  (s)     = 
u  k^    '  2, 


S    (S  -  Sj^) 

By  a  procedure  completely  analogous  to  that  just  described  for  the  radial  deflec- 


tion, the  results  can  be  written 

i(ps 


f 


— ds     ^     7riRes(0)  +27riRes[    f,  (s,  )] ,  0>O 

^  (^  -  \)  b^  >  0 

k 


(5-33) 
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/-     ei(0+d)s 
J     ~ -ds     =     7riRes'(0)  +  27riRes[^f|^(s^)l,  (p^  >  0         (5-34) 


S    (S   -   Sj^) 


b,    >  0 
k 


where 


-\(l>  +  ^\4> 


^^«fufk(\)]      = r  (5-35) 


-bj^(0+d)+iaj^(0-+d) 

^^«Iu^k(\)J      ^ i (5-36) 

i\  +  i  bj^) 

To  find  the  residue  of  the  functions  at  s  =  0  they  are  expanded  in  a  Laurent  series 
about  the  point  s  =  0 . 


e 


=     J-  +  ii  _  ^   _  i*^    s    ^ 


s2  g2         s  2!  3! 


-1-     =     -  J-     1   .  ^   .  ^   , 


e""~  _  1   /    1     ,  \  ^  -•  .         .2 


i0s  /  i<^  ^ 


s^(s-sj^)  \\s^  s  s/         \         2 

The  residues  can  be  evaluated  as 

Res(O)     =     -  -1-   _  i  A  ^3_3^j 

\  k 


Res'(O)    =     -  J-   _  i  ±±A 


\'  '^ 


(5-38) 
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so  that 


/ 


i(j)S 


s    (s  -  Sj^) 


as.    -f^(t-*)--5 


b,    >  0 
k 

0  >  0 


/ 


i((/)+d)s 


ds 


s    (s  -  Sj^) 


-'tk'^H 


i((/)^)S, 


+  i((p+d)      +  27ri 


.   e 


b,    >  0 
k 


k  d  +(p  >  0 

Integrating  in  the  lower  half  plane  for  <^  <  0,  Equation  (5-30)  can  now  be  written 
in  the  form 

i<^s. 


F,  (0)     =     i  -^-^ 
u    k^^'  s, 


-nt-*)^^ 


-i(t^'<*^^')^^ 


H(0) 
i(0-hd)s,^ 


H(0+d) 


+  — 


,/l+i0s^ 


2  \        s, 


H(-0) 


^  / 1  +  i(cp-Hi)  s^ 


H(-0-d)  1  ,  b,    >  0  (5-39) 


Performing  the  integrations  for  b,    <  0,  the  final  result  can  be  written  as 


F,  (0) 


sgn(bj^)ic^aj^ 


[-2(l+i<^sj^)-fe       ^J 


X  H[sgn(bj^)0]   +  |(1  +i  0  Sj^)  H[-sgn(bj^)0] 

,  i((/)-Kl)s    1 

■■^[l+i{(p-Hi)s^]+e  '^    H[sgn(bj^)(0-kl)] 


+  i-  [l+i(0+d)Sj^]  H[-sgn(bj^)(0+d)] 


(5-40) 
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The  solutions  given  by  Equations  (5-7),  (5-19),  and  (5-31)  with  the  correspond- 
ing F.    functions  (5-13),  (5-27),  and  (5-40)  are  substituted  into  the  governing 
equations  of  motion  in  Appendix  G  to  show  that  the  solutions  satisfy  the  differen- 
tial equations. 

Summary  of  Deflection  Expressions 


The  following  expressions  summarize  the  solution  for  the  deflections. 

4 
2    Fj^(0)  (5-41a) 


k=l 
4 


m^ 


k=l 


(5-41b) 


}M1    = 


.^^    ±i^   +Y       F,(0) 
c         q  Z^    u    k^^' 


(5-410) 


9  'O 


k=l 


where 


F^ict» 


F,  (0) 


iSj^(0+d) 


(     IS     <p  1 

sgn(bj^)iCgaj^     e     *"      H[sgn(bj^)0]  -  e     '^  H[sgn(bj^)(0+d)] 

(5-41d) 


sgn(bj^)  /3j^ 


1   ^    '%"> 
-2   +e 


H[sgn(bj^)0]   +2  H[-sgn(bj^)0] 


^        is,((^+d)' 


H[sgn(bj^)((^+d)] 


+ 1  H[-sgn(bj^)(<^+tl)] 


(5-41e) 
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F,  (0) 
u    k^^ 


sgn(bj^)ic^aj^ 


[-1 


(1  +is^0)  +e 


i  ct>  Sj^-i 


H[sgn(bj^)0] 


+  i-(l  +iSj^</,)H[-sgn(bj^)0] 


-|  [1  +is^(0-Hl)]    +e     "^  H[sgn(b^)(04d)] 


2    [1  +isj^(<^+d)]   H[-sgn(bj^)(0+d)] 


(5-41f) 


sgn(bj^)     = 


1 ,  bj^  >  0 


-1,  bj^<0 


H(y) 


y  <  0 

y  >i 


a 


Px 


c     s,     +  c 

7      k  8 


k  4  '  ^k  4 

n  (\-v  n  <^k-v 


m=l 


m=l 


kT^m 


s,  =  a.    +  i  b,   (k  =  1,  , . . ,  4)  are  the  roots  of  the  characteristic  equation 


D(s)     =     c     s    +ic     s^+c    s^+ic     s+c       =     0 

4  3  2  1  O 


57 
and  the  coefficients  are  defined  as 


c^     =     -exe  f 

1  3  8 


C 

2 


=     f  (e  c  -  e  e  )  +e  f  (f    -  rx^) 

7      O  12  3   6     2 


C  -      -exf  C 

3  7 

c^  =     (f^  -  rx^)f  c 

4  3  7 


c'      =     f  (e    +  e  ) 

5  7^     3  l' 


c^      =     e  e 

6  13 


<    =    fA^V^) 


c'      =     e  f 

e  3  e 


c'      =     e  f    +c' 

9  17  8 


e=E+Pr  f=l  +  F 

o              o          o  o 

e       -     E     -  P  f,     =     1  +  F  -  D 

■L                       -L  1                                               1 

e       =     E     -  P  +  2(G  +  Pr^)  f       =     F  -  Pr 

2                    1  0                       2                                  o 

e=G  +  Pr  f=G  +  F 

J                                  O  3 

e.     =     I    -  r,  f       =     D  r     +  I  f 

^Ol  4                     1100 

®5     "     ^n^  ■  ®.^^  f      =     r  f    +rf 

5                    0112  5                    120 


(5-42) 
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The  derivatives  of  the  deflections  are 


4 
k=l 

% 

4 

1_      1 

q  c  q  ^-i^    u    k^^  ^  ' 


c 

■» 

"9  "O 


k=l 


where 

^k 

f;(0)  = ^—  F,  (0)  = 

9     k 
u^k<^)      =     -^w^k(^)     = V-,.F,«^)  (5-44) 


2     , 

c     s,      +  c 

7     k           e 

c 

9 

Fk(^) 

C 

6 

w\«^) 

2 
C      S,      +C 

7      k            8 

</>,  <^+d 

/ 

0 

Solution  for  No  External  Damping 

For  the  case  of  no  external  damping,  the  solutions  as  given  in  Equation  (5-41) 
are  not  directly  applicable.    It  will  be  noticed  that  the  forms  of  the  solutions  are 
dependent  upon  the  signs  of  the  imaginary  parts  of  the  complex  roots.    Figures 
4 . 2  through  4 . 8  show  that  there  are  regions  in  which  there  are  only  real  roots 
of  the  Characteristic  Equation  with  no  damping,  thus  causing  a  problem  of  non- 
uniqueness  of  the  solutions . 

Following  the  method  of  Achenbach  and  Sun  (13)  the  undamped  solution  will 
be  obtained  uniquely  by  assuming  the  undamped  solution  is  the  limit  of  the  damped 
solution  as  the  damping  approaches  zero.    In  this  manner  the  sgn  functions  in  the 


59 

F^  functions  for  Equation  (5-41)  can  be  determined.     Figure  5.4  provides  an 
example  of  the  behavior  of  a  set  of  roots   as  the  damping,    e,    approaches 
zero.    This  figure  shows  the  type  of  the  four  roots  (s..  =  root  i  with  damp- 
ing €.)  for  heavy  damping  to  be  two  complex  and  two  imaginary,    root  1 
being  very  near  the  origin.    As  the  damping  is  decreased,  the  imaginary  roots 
approach  one  another  and  finally  meet  and  separate  which  gives  four  complex 
roots.    Meanwhile  the  other  complex  roots  are  also  approaching  the  real  axis. 
This  establishes  the  correct  sign  for  the  imaginary  part  of  each  root  in  the  limit 
as  c  -*  0  and  the  roots  all  approach  the  real  axis. 

Form  of  the  Radial  Deflection  in  Region  IV 

The  roots  are  all  real,  having  the  form  ±a   ,  ±a   .    From  Equation  (5-41b), 
the  radial  deflection  expression  (when  the  proper  signs  are  established)   in 
Region  IV  becomes 

O  L       1  1  2  2  J 

/3     r  ia    0  la  {<p+d)  1 

-  ^     e     ^      H(-0)  -  e     ^  H(-cp-d) 

(3     r  -ia    0  -ia  (04d)  "l 

+  —     e  H(-0)  -e  H(-0-d) 

[ia    </)  ia  {(p+d)  "I 

e     ^     H(<^)  -  e     ^  H(<p+d) 

p     r  -ia.    (p  -ia  {(p+d)  1 

-^     e       ^      H(0)-e       ^  H((^-H1)  (5-45) 


a 

2 
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23 


s. .  =  Root  i  for 
External  Damping 

€ 

e 
o 

=      0 

^1 

=      0.1 

^2 

=       1.0 

S 

=    10.0 

Figure  5.4.    Loci  of  the  Roots  of  the  Characteristic  Equation 
as  the  Damping  Approaches  Zero 


61 


Expanding  the  /?.  coefficients  gives 


/?. 


-/?, 


c    a    +  c 

7     1         e 

2a  (a^  -  a^) 

1       1  2 


c    a    +  c 

'^^-     2         e 

2a  (a^  -  a") 

2'     1  2' 


Substituting  expressions  (5-46)  into  (5-45)  gives 


(5-46) 


W 

q 


o 


ca+c       l/ia(;&        -ia^ 

2a^(a^  -  aT) 
11        2 


[ia  (0+d)        -ia  (0+d)' 
e     ^  +e       ^ 


H(-0-d) 


c     a    +c       l/iad)        -iad)\ 
7     2        e     1/^     2^^^       ^         H(0) 


22        ? 

2a  (a    -  a  ) 

2'   1        2' 


•[• 


ia^(0+d)        -ia^(0+d) 


+  e 


] 


H(0+d) 


o2  „2 

a    a 

1    2 


[H(0^)  -  H(0)] 


or 


2  2 

a    a 

1    2 


[H(0-Hi)  -  H(0)]    - 


ca+c 

7     1        e 

a^(a^  -  a^) 

1^  1        2' 


^  ["cos  (a^0)  H(-0)  -  cos  [a  (0+d)]  H(-0-d)l 


c    a^  +  c 

7    2 


af (a^  -  a") 

2^    1  2' 


cos  (a  0)  H(0)  -  cos  [a  (0+d)]  H(0+d)|  (5-47) 

-  a^\  L  2  2  J 
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Therefore,  for  the  three  distinct  regions  of  the  cylinder,  the  solutions  are: 
Solution  behind  the  load  (0  <  -d) 

f"  =     -  ^"7 f  {cos  (ad)) -cos  [a  (0+d)]}  (5-48) 

1^  1        2' 

Solution  under  load  (-d  <  0  <  0) 

XT,  c  ca+c  ca+c 

--  =     -  — ^- ^— i cos  (a  0)  +      ^    ^ 2.  cos  [a  ((/)+d)]       (5-49) 

q                  22  2,  2        2,          ^   1^'         2,  2        2,          '  2^'^     "       ^         ' 

o             a    a  a  (a    -  a  )                        a  (a    -  a  ) 

12  i''  1       2'                          2^  1        2' 

Solution  ahead  of  load  (0  >  0) 

[cos  (a^0)  -  cos  [a^(0+d)  ]|  (5-50) 


2  ^ 

Wc    a    +  c 
_  7       2  8 

%  a^(a^  -  a^) 

2^  1        2' 


Form  of  the  Radial  Deflection  in  Region  VII 

As  in  Region  IV,  the  roots  are  ±a   ,  ±a   .    However,  there  can  be  no  deflection 
ahead  of  the  load  in  this  region  because  the  load  speed  is  greater  than  any  of  the 
wave  speeds  in  the  material.    This  zero  displacement  comes  about  mathematically 
because  the  roots  all  approach  the  real  axis  from  the  negative  imaginary  direction. 
The  solution  for  this  region  has  the  form 

TTT  0  ca+c 

W  e  r^^,  J,  TX,  J  1X-,  T  1  8 


[H(-0)  -  H(-0-d)]    - 


q                   22^           \    T      /J  2,  2        2, 

o              a    a  a  (a    -  a  ) 

12  1^  1        2' 

X  /cos  (a  0)H(-0)  -  cos  [a  {(p+d)  ]H(-0-d)l 


2  , 
ca+c 

+         7       2  8 

2      2  2, 

a  (a    -  a  ) 
2^  1        2 


icos  (a^0)H(-0)  -  cos  [a^(0+d)]H(=0-d)l  (5-51) 
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Form  of  the  Radial  Deflection  in  Region  V 

The  roots  in  this  region  have  the  form  ±a,  ±ib.    The  radial  deflection  in  this 
region  is  given  by  the  expression 


W  ^e 


%         a^  b^ 


[H(-0)  -  H(-0-d)] 


c    a    +  c 

7 
2      2 

a  (a 


1-  fcos(a0)H(-0)-cos  [a(^+d)]  H(-(^-d)  1 

+  b  )  L  -J 


u2     . 

c    b    +  c 


^      ;    a       I     fe"''*^  H(c^)  -  e^^  H(-<^)  -  e'^^^^^)  H(0^) 
2b  (a    +  b  )    •- 

+  e^('^"^)H(-0-d)]  (5-52) 

Form  of  the  Radial  Deflection  in  Region  VI 

In  this  region  the  roots  are  all  imaginary,  of  the  form  ±ib   ,  ±ib   .    This 
gives  an  exponentially  decaying  solution  as  in  Regions  I  and  III.    The  radial  de- 
flection expression  is  given  below. 


e 


^    =      -  -T^    [H((/)-Hl)  -  H(0)] 
^o  b    b 

1      2 


c     -  c    b^      r  -b    0  b    0  -b  {(p+d) 

_!_     e        -     ^  ^ 


^-^  \e     ^      H(0)  -  e   ^      H(-0)  -  e     ^  H{(p+d) 

+  6   ^  H(-0-d)J 

r   b„  0  -b    <^  u  vv"; 

e  ^      H(-0)  -  e  H(0)  -  e  H{-<p-d) 

J<P^)  -1 

+  e  H(04td)  (5-53) 


2b^(b^  -  b  ) 

^    ^        ^  b  (<^^) 


Cg-c^b^     r  b^0  -b,<^  b^(0+d) 


2bp(bi  -  b^) 

-b^((^+d) 
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Comparison  of  Solution  with  Other  Results 

As  a  comparison  of  the  results  of  this  analysis  with  another  theory,  the 
static  problem  of  a  distributed  pressure  load  on  an  isotropic  shell  was  considered 
as  shown  in  Figure  5.5, 


-• i' 


T 

R 


K^  ^2  ^ 


Figure  5 . 5    Static  Load  Problem 

For  the  static  problem  shown  above  the  roots  of  the  characteristic  equation, 
excluding  prestress,  are  complex.    The  solution  for  the  region  under  the  load 
as  given  by  the  present  theory  can  be  reduced  to 

—     =     K       2  -  e  cos  aL    -  e  cos  aL 


-a  L  -a  L 

+  K    I  e  sin  a  L    +  e  sin  a  L 

2   \  2  1 


(5-54) 


where 


^a     =     -f- 


\     -    t 


65 
This  problem  is  solved  by  Timoshenko  (14)  and  the  deflection  given  by  his 
theory,  when  put  in  a  compatible  form,  becomes 

W  (  -^^2  -^^1 


—     =     KM2-e  cos /3i^  -  e  cos  (3£^j  (5-55) 

The  second  term  appearing  in  Equation  (5-54)  is  missing  from  Equation 
(5-55).    This  additional  term  arises  because  of  the  inclusion  of  shear  deflection 
which  was  not  present  in  the  Timoshenko  theory,    A  comparison  of  the  results 
of  these  deflection  expressions  is  made  in  Figure  5.7. 

Numerical  Results 

A  computer  program  was  developed  for  the  calculation  of  the  displacements 
and  stresses  determined  in  this  research.    The  general  expressions  given  by 
Equations  (5-41)  were  programmed  for  the  displacement  solutions  and  the  stress 
calculations  are  discussed  in  Chapter  VI.    A  flow  diagram  of  the  computer  pro- 
gram is  shown  in  Figure  5.6.    The  details  of  the  program  can  be  found  in  Ap- 
pendix H.    It  is  written  in  Fortran  for  time  share  computer  application. 

Comparison  of  Results  with  Other  Solutions  for  a  Static  Load 

The  radial  deflection  for  a  static  distributed  load  on  a  cylindrical  shell  is 
given  in  Reference  (14),    As  a  check  on  the  solution  this  static  problem  was 
solved  using  the  present  results  and  the  comparison  is  shown  in  Figure  5.7, 
The  results  agree  very  well.    The  effect  of  variation  in  the  thickness-to- radius 
is  also  illustrated  in  Figure  5.7.    The  rotation  and  axial  deflection  are  shown 
for  this  static  problem  in  Figure  5.8.    In  addition  to  showing  the  form  of  the 
displacements  for  the  static  load,  Figures  5.7  and  5.8  serve  as  a  basis  against 
which  the  dynamic  displacements  can  be  compared.    The  deflections  are  sym- 
metric about  (p  =  -0.5  for  the  static  load. 
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Figure  5.6.     Flow  Diagram  for  Computer  Program  for  Deflection 
and  Stress  Calculations 
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Figure  5.7.    Radial  Deflection  for  Static  Load  on  an  Isotropic  Shell  (fi  =  0.3) 
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Figure  5.8.    Displacements  for  a  Static  Load  on  an  Isotropic  Shell 
(/i=  0.3,  h/R  =  0.01) 


69 
Summary  of  Deflection  Response  for  Shells  under  Various  Load  Velocities 

A  summary  of  some  of  the  types  of  deflection  patterns  assumed  by  a  shell 
for  increasing  load  speed  is  shown  in  Figure  5.9.    For  the  particular  properties 
used  for  this  example,  the  various  regions  (root  types)  associated  with  each 
waveform  can  be  found  by  inspection  of  Figures  4.2  through  4.8.    For  example, 
for  no  damping,  positive  prestress  corresponding  to  internal  pressure, 
h/R  =  0,001  and  the  material  properties  given  in  Figure  5.9  (properties  are 
those  corresponding  to  an  isotropic  shell  as  shown  in  Appendix  I),  Figure  4.2 
can  be  used  to  associate  load  speed  with  root  type. 

Following  the  vertical  line  of  h/R  =  0.001,  it  is  evident  that  A    =  1  lies  in 
Region  I  where  the  roots  are  all  imaginary.    This  gives  a  critically  damped 
exponentially  decaying  solution  as  shown  by  Equation  (5-53)  and  is  shown  in 
Figure  5.9(a).    As  the  load  speed  increases  Region  II  is  entered  where  the  roots 
are  complex.    This  is  the  form  of  the  static  load  problem  roots,  and  if  no  pre- 
stress existed  Region  n  would  extend  from  zero  load  speed  up  to  the  first  criti- 
cal,  which  is  at  A    .    The  solution  for  X    =  2  is  shown  in  Figure  5.9(b),  and  is 
exponentially  decaying. 

The  response  becomes  sinusoidal  after  crossing  \  ~  X    .    At  a  load  speed 
just  greater  than  A    the  deflection  response  has  a  very  short  period.    A  small 
amplitude  wave  train  precedes  the  load  and  a  large  amplitude  wave  follows  it. 
As  the  load  speed  increases  the  sine  wave  period  increases  as  shown  in  Fig- 

2 

ure  5.9(d)  for  A    =30.    These  sinusoidal  deflection  patterns  are  in  Region  IV 
where  the  roots  of  the  Characteristic  Equation  are  all  real .    The  mathematical 
expression  forW/q    is  given  by  Equation  (5-47).    Crossing  A    into  Region  V,  the 
roots  are  real  and  imaginary.    Equation  (5-52)  gives  the  radial  deflection,  and 

2 

Figure  5.9(e)  shows  the  response  to  be  a  long  period  sine  function  for  A    =  500. 
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Figure  5.9.    Radial  Deflection  Shape  for  Various  Load  Velocities 
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Jumping  to  Region  VII  brings  a  longer  period  sinusoidal  oscillation  as  illus- 
trated by  Figure  5.9(f)  and  the  response  in  this  region,  where  the  roots  are  again 
all  real,  was  discussed  previously.  The  radial  deflection  is  given  by  Equation  (5-51). 

Because  Region  VI  covers  such  a  limited  range  in  velocity  the  response  was 
not  included  in  the  summary  but  is  discussed  later. 

Region  II  Hcsi)onse 

A  study  of  the  response  of  an  isotropic  shell  at  a  load  speed  below  the  first 
critical  was  made  to  determine  the  effect  of  external  damping.    These  results 
are  shown  in  Figure  5. 10  where  the  damping  ranges  from  very  light  to  very  heavy. 
Of  course,  when  damping  is  introduced  the  root  form  is  no  longer  the  same  as 
that  of  Region  II. 

Region  IV  Response 

2 

The  short  period  sinusoidal  response  of  the  radial  deflection  at  A    =  2 . 7  is 
shown  in  Figure  5.11.    As  the  radial  damping  is  increased  this  response  is 
changed  drastically  as  shown  in  Figure  5.12.    The  response  for  a  damped  sys- 
tem, which  was  in  Region  IV  with  e  =  0,  approaches  closely  that  of  the  Region  II 

2 

behavior.    Figure  5.13  shows  the  radial  response  at  X    =5  and  10.    The  maxi- 
mum amplitude  remains  constant  as  the  period  of  the  wave  increases  for  greater 
load  velocities. 

Damping  Effect  on  Regions  V  and  VU  Response 


2 

The  effect  of  damping  on  the  wave  forms  for  X    =  500  and  2000  is  shown  in 
Figures  5.14  and  5.15.    The  amplitudes  of  the  sinusoidal  deflection  response  are 
initially  decreased,  and,  as  the  damping  becomes  greater,  the  response  becomes 
critically  damped  and  the  deflection  approaches  zero  with  an  increase  in  distance 
from  the  load. 
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Figure  5. 10.    Radial  Deflection  Response  for  Variations  in 
Radial  Damping  (x2  =  2,0) 
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Figure  5.11.    Radial  Deflection  Pattern  Immediately  Above  the  First  Critical 
Load  Speed 
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Figure  5. 12.    Change  in  the  Radial  Deflection  Pattern  with  Increasing 
Damping  (A^  =  2 . 7) 
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Figure  5. 13.    Deflection  Wave  Form  at  X^  =  5  and  10 
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Figure  5. 15.    Effect  of  Damping  for  X^  =  2000 
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Region  VI  Response 

This  region  has  a  deflection  pattern  which  is  almost  totally  axial.    Three 
critical  velocities  have  been  crossed  to  get  into  this  velocity  range  which  corre- 
spond to  the  longitudinal ,  shear ,  and  bar  wave  speeds  .    Therefore ,  there  can 
be  no  bending  effect  transmitted.    The  behavior  is  like  an  axial  compression  on 
a  membrane  which  expands  radially,  as  shown  by  Figure  5 .  16 .    The  effect  of 
prestress  is  observable  in  this  figure.     The  maximum  deflection  is  increased 
by  about  20  percent  when  going  from  external  hydrostatic  pressure  to  internal 
hydrostatic  pressure. 

This  axial  mode  of  deflection  also  appears  in  other  velocity  ranges.    For 
instance,  the  broken  lines  in  Figure  5.16  show  the  behavior  at  \    =  1001.    This 
is  the  range  between  A    =  1000  and  X    =  1002  where  the  roots  are  real  and 
imaginary.    Another  example  is  shown  in  Figure  5.22  where  En  /E^    is  less 
than  0.08  and  this  occurs  in  Region  VUI  as  shown  in  Figure  4.4, 


2 

Deflection  Behavior  in  the  Vicinity  of  An 


CR 


2  _   ,2 


The  first  (lowest)  critical  load  speed  occurs  at  Aip,„  =  A    =  2.552.    Figure 
5.17  illustrates  the  unbounded  response  of  the  deflection  as  that  speed  is  ap- 
proached .    The  effect  of  damping  on  the  maximum  deflection  is  also  illustrated . 
There  are  four  other  critical  velocities  as  discussed  in  Chapter  IV. 

2 

Effect  of  Prestress  on  A.  ^r^ 

^^CK 

The  effect  of  the  axial  prestress  on  the  location  of  the  first  critical  load 
speed  is  shown  in  Figure  5.18.    This  effect  is  also  observable  in  Figure  4.8, 
since  the  first  critical  load  speed  is  at  A^ .    The  circumferential  prestress  does 
not  have  a  significant  effect  on  load  speed,  as  shown  in  Figure  4.7. 


79 


_u 


Figure  5.  16.    Deflection  Response  for  Variations  in  Circumferential 
and  Axial  Prestress 
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Superposition  of  Step  Loads 

The  effect  on  the  radial  deflection  can  be  observed  in  Figure  5 .  19(a)-(f) 
where  the  load  length  was  varied  from  0.1  to  5.0.     By  superposing  various  com- 
binations of  step  loads  it  is  possible  to  approximate  any  shape  of  load  desired. 
As  an  example  of  this  type  of  application,  the  radial  deflection  response  from 
a  symmetric  sine  wave  type  load  and  a  sharp  edged  pressure  front  was  deter- 
mined.   The  results  of  these  calculations  are  presented  in  Figures  5.20  and 
5.21,  respectively. 

Study  of  Material  Properties  Variations 

A  look  at  the  effect  of  decreasing  the  Eq  /E^    ratio  is  summarized  in  Fig- 
ure 5.22.    Starting  in  Region  IV,  as  can  be  seen  in  Figure  4.4,  the  ratio  is  de- 
creased from  1.0  (as  for  an  isotropic  material)  to  0.04.    As  the  ratio  is  lowered, 
the  maximum  deflection  gets  large  rapidly,  and  becomes  unbounded  as  ?^^  =^  ^IpR 
is  approached.    After  crossing  ^^cn  ^"^°  Region  VIII,  the  strength  in  the  cir- 
cumferential direction  is  of  course  very  low  and  the  material  is  of  little  interest 
for  engineering  applications. 

The  same  type  of  response  as  in  Figure  5.22  will  be  obtained  by  increasing 
E,,  /Ej.    significantly.    This  can  be  observed  by  inspecting  Figure  4.6.    For 
other  types  of  material  property  variations,  the  general  response  can  be  pin- 
pointed by  observing  the  type  of  roots  at  the  particular  location  through  the  use 
of  Figures  4 . 2  through  4 . 8  and  using  the  numerical  results  presented  here  show- 
ing similar  calculations  of  deflections . 
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Figure  5.22.     Radial  Deflection  Response  for  Variations  in  the 
Circumferential  Modulus 


CHAPTER  VI 
STRESSES 

Development  of  Stress  Equations 


For  the  axisymmetric  loading  the  stress  and  moment  resultants  become 

N       =E^+^M   +  E^  (6-1) 

X  X  ax       R    ax        I'  R  ^     ' 


N=(^E+^^^   +   e|^  (6-2) 


^x  =  ^  i  ^  °x  i  <^-=' 


«e'-    -  X  I  ^  '^^  If  ('^-^) 


Assuming  the  bending  stress  is  linear  across  the  thickness  of  the  shell,  the 
stresses  can  be  put  in  the  nondimensional  form 


(6-6) 


E  d<^  1  d4>       ^1^  E^ 

X  ^  ^  ^o 


N^  JTT  O'fl 

^     =      E    W  +   E     ^     =     ^  (6-7) 

E  o  1  dd)  Ev 
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"^x    _    du  ^  di^         o  R  !^ 

X  ^  ^  Xq 

GQ 


X 


G  ^V^       d(p 


The  stresses  are  therefore 


„  /  ,   ^  dW  \  ^xz 


HL.  =  ^  ym  ^^  Him  ^uim  (6_n) 


^Xo%  ^        %  '        %  ^ 


o 


^xo  %  °      %  ^       %  ^        ^ 


-!^^    =    i  Ji/UM  ^i!!^ 


Ey    q  2    R  I     q  q 


Ev  q 

^o  ^o 

2    R  \     1      q,       - 

Ev  q 

Xo  ^0 

(6-1.3) 


(6-14) 


(6-15) 
The  stresses  at  the  outer  and  inner  surfaces  of  the  shell  are  now  determined. 

(6-16) 


X  outer 
q    E-, 

^O        Xq 

q  Ex 

^0         Xq 

+ 

q  Ex 

^O         Xq 

^""o^outer 

q   E„ 

^o     ^o 

^0m 

q  Ex 

^o       Xq 

+ 

""eh 
q   Ex 

^o     ^o 

(6-17) 
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(a  ). 
X  inner 

q     Ex 

^^b 
q     Ex 

^^e'inner 

q  Ex 

"^m 
%% 

q  Ex 

(6-18) 


(6-19) 

Both  the  individual  stresses  given  by  Equations  (6-11)  -  (6-15)  and  (6-16  - 
(6-19)  are  calculated  by  the  computer  program  which  is  outlined  in  Figure  5.6. 
The  program  listing  is  given  in  Appendix  H. 

There  are  various  relationships  which  exist  between  parameters  used  in 
this  research.    Some  of  these  relationships  are  collected  in  Appendix  I  for  refer- 
ence purposes,  for  use  in  obtaining  the  dimensional  form  of  stresses  and 
deflections . 

Numerical  Results 

For  an  isotropic  shell  under  a  static  ring  load ,  P ,  the  second  derivative  of 
the  radial  deflection  is  given  in  Reference  (14)  as 

d     w  2  P 

=     2/3     e  (sin  /3x  -  cos  /3x)  (6-20) 

dx^  8/3^  D 

The  solution  for  the  static  problem  shown  in  Figure  5.5  can  be  obtained  by  re- 
placing P  by  q  dx  and  integrating  so  that 


dx  *J|  8/3     D 


l^     =    J      .«=  _3^   /^^ 


2  /  ^^        ~3 ^  ^^^"^   ^'^   ~   ^°^    l^^^ 


K 


C  2     a  dx       '^^ 

+   I        2/3        ^  e  (sin  ^Sx  -  cos  )3x)  (6-21) 

•J  8/3=  D 
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The  moment,  M    ,  in  Reference  (14)  is  assumed  to  be 

,2 

M       =     -  D  ^-f  (6-22) 

Combining  Equations  (6-21)  and  (6-22)  and  assuming  a  linear  stress  distri- 
bution through  the  thickness,  the  bending  stress  is  determined. 

<^Xb  3  /    -^^2  -^^1  \ 

'e  sin(S£     +e  sin  /?  £    I  (6-23) 


2^3(1  -  M^) 


where 


4     ^      3(1  -  M^) 

2       2 

R    h 


(6-24) 


The  stress  given  by  Equation  (6-23)  is  shown  in  Figure  6.1  in  comparison 
with  the  results  derived  in  this  work.    The  results  given  by  (6-23)  are  good  for 
thin  shells  and  get  worse  as  the  shell  becomes  thicker.    The  axial  and  tangential 
surface  stresses  are  shown  in  Figure  6.2  for  the  static  load  problem.    The  com- 
puter program  developed  can  be  used  to  calculate  the  individual  stresses  and/or 
the  surface  stresses  as  desired. 
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0         (p      x/R 


R 


Present  Theory 
Equation  (6-23) 


x/R 


Figure  6.1.    Bending  Stress  in  an  Isotropic  Shell  Under 
a  Static  Load  (/i  =  0.3,  d  =  1) 
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<t>      x/H 


Figure  6.2.    Surface  Stresses  in  an  Isotropic  Shell  Under 
a  Static  Load  (;i  =  0.3.  h/R  =  0.1) 


CHAPTER  VII 
CONCLUDING  REMARKS 

Conclusions 

An  analysis  has  been  made  of  an  orthotropic  infinite  circular  cylinder  sub- 
jected to  an  axisymmetric  discontinuous  pressure  loading  of  finite  length  moving 
in  the  axial  direction  with  a  constant  velocity.     Fourier  transforms  were  used  to 
solve  for  the  steady  state  solution  and  proved  to  be  very  useful  for  this  type  of 
application.    The  problem  of  nonuniqueness-of-solution  arose  for  the  case  of  an 
undamped  cylinder .    This  was  handled  by  taking  the  undamped  solution  as  the 
limit  of  the  damped  solution  when  the  damping  approached  zero . 

Investigation  of  the  higher  load  speeds  showed  that  there  are  five  critical 
load  speeds,  the  two  highest  being  very  close  together  and  corresponding  approxi- 
mately to  the  plate  wave  speed.    Lowering  of  the  tangential  modulus  ratio  was 
found  to  cause  a  rapid  increase  in  deflections .    An  increase  of  speed  above  the 
bar  wave  speed  was  found  to  cause  a  mode  of  deflection  which  is  mainly  axial 
in  nature  as  opposed  to  the  low  speed  bending  type  deflection  pattern.    The  method 
of  superposition  was  used  to  demonstrate  the  adaptability  of  the  results  to  other 
forms  of  loads. 

Although  the  shell  theory  used  was  linear ,  small  deflection  theory ,  the 
results  are  useful  for  loads  in  the  range  for  which  most  structures  are  designed. 
The  general  response  picture  can  be  obtained  from  the  results  presented  here 
and  the  location  of  the  critical  load  speeds  observed .    The  computer  program 
developed  here  can  be  used  to  investigate  in  detail  particular  problems  of  in- 
terest.   Because  of  the  increasing  frequency  of  applications  of  cylindrical  shells 
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in  dynamic  loading  conditions  in  military,  government,   and  industrial  programs, 
the  critical  load  speed  environment  should  be  considered  in  the  design  require- 
ments.   Programs  utilizing  thin  walled  launching  tubes  for  launching  missiles  are 
a  specific  example  of  such  applications. 

Suggestions  for  Future  Work 

For  large  deflections  the  theory  used  in  this  analysis  becomes  invalid.    A 
natural  extension  of  this  work  would  be  the  investigation  of  the  stability  through 
use  of  large  deflection  theory.    Also  the  increased  use  of  composite  materials 
for  structural  applications  makes  the  inclusion  of  anisotropic  material  a  worth- 
while endeavor .    In  this  analysis  the  axial  and  circumferential  prestress  were 
considered.    Another  parameter  which  may  prove  to  have  a  significant  effect  is 
that  of  torsional  prestress. 

The  computer  program  developed  for  the  calculation  of  stresses  and  deflec- 
tions in  the  shell  is  by  no  means  efficiently  programmed.    Extensive  use  would 
warrant  an  investigation  to  reduce  run  time. 


APPENDIX  A 
FOURIER  TRANSFORM  OF  THE  FORCING  FUNCTION 


Let  L  denote  the  transform  of  the  pressure  function  given  on  the  right  hand 
side  of  Equation  (2-34b);   the  integral  to  be  evaluated  is 


^f     =      f    % 


[H(-R0)  -  H(-R0  -  Vt  )]   e"^^*^  d0 


(A-1) 


or 


%/ 


[H(-</))  -  H(-</)  -  d)]   e"^^*^  d0 


(A-2) 


where 


d     = 


Vt 
o 

R 


{A-3) 


H(-0) 


i: 


<;&  >  0 
0  <  0 


(A-4) 


H(-(/)-d)     - 


c 


0  +  d  >  0 


(/)  +  d  <  0 

Integral  (A-2)  therefore  can  be  broken  up  into  two  integrals, 
0  -d 


I 


f 


f  e-'^^d4>  -    f     e-' 


-is* 


d0 


(A-5) 


(A-6) 
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Furthermore ,  the  second  integral  can  be  written 

-d  0  0 

f     e-^^'^dc/)    =     f    e'^'^d<p  -   f   e-^"<^d0  (A-7) 

—  00  _00  — (j 

Therefore, 


I,    =     q 


„/e-'-a*=^[e-'-];^ 


T      _     _o    ,1        isd.  ^o   ,  isd      ,.  ,.    „, 

If     -     Iil(l-e       )     =     -^(e         -1)  (A-8) 

Equation  (A-8)  gives  the  transform  of  the  pressure  loading. 


APPENDIX  B 
SOLUTION  OF  EQUATIONS  FOR  THE  TRANSFORMED  DEFLECTIONS 

From  Equation  (3-7)  the  set  of  equations  to  be  solved  is 


(l+F-r^X^)s^ 

-i(E^-P)s 

(D^-I„A^)s^ 

U(s) 

i(E^-P)s 

(G+F-rA^)s^-ie\s 

+E  +Pr 

o         o 

-i(G+Pr^)s 

W(s) 

(D  -I  X^)s^ 

^    1     o      ' 

i(G+Pr  )s 

+G+Pr^   _ 

m 

^o  ,       isd 


(B-l) 
First,  the  determinant  of  the  coefficients,  denoted  as  A,  will  be  determined. 


A     = 


1+F-r  X 

1 


i(E^-P)s 


(D  -I  A^)s^ 

^    1     o      ' 


■i(E  -P)s 


(G+F-rA^)s^-ieAs+E  +Pr 

^  '  GO 


i(G+Pr  )s 

o 


(D  -I  A^)s^ 
^    1     o      ' 


■i(G+Pr^)s 


(D  -I  A^)s^+G+Pr 

^    1     o      '  o 


(B-2) 


Expanding  this  determinant  gives 


2      2 


2       4 


2        3 


A     =     (1  +  F  -  r^A  )s    ^  (G  +  F  -  rA  )  (D^  -  I^A  )s     -  ieA(D^  -  I^A  ) 


(D^-I^A  )(E^+Pr^)+(G+F-rA  )(G+Pr^) 


-  ieA(G  +  Pr^)s  +  (E^+  Pr^)  (G  +  Pr^) 


2.      4 


2.2 


2.     6 


2     4 


-  2(E^  -  P)  (G  +  Pr^)  (D^  -  I^A  )s     -  (D^  -  I^A  )    (G  +  F  -  rA  )s 
+  ieA(D^  -I^aYs^-  (Eq  +Pr^)  (D^  -I^A^s^  "  (1  +  F-  r^A^)  (G+Pr^)"s 

-  (D^  -  I^A^)  (E^  -  P)^s*  -  (G  +  Pr^)  (E^  -  P)^s^  (B-3) 
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Collecting  coefficients  on  powers  of  s  gives  the  expression  for  the  deter- 
minant in  the  form 


A     =     (G  +  F-  rA^)  (D^  -I^A^)  ^(1  +  F-  r^\^)  -(D^  -I^\^) 
+  iekCD^  -  l^X^)  I  (D^  -  I^\^)  -  (1  +  F  -  r^A^)l  s^ 


(1  +  F  -  r  A" )  -  (D    -I  X^) 

1      '       ^    1        o     ' 


+  i  (E    +  Pr  )  (D    -  I  A  ) 

'  ^    o  o'  ^    1        o     ' 


+  (1  +  F  -  r^\  )  (G  +  Pr^) 


(G  +  F  -  rA  )  -  (G  +  Pr  ) 


-(D^-Io^*^) 


(E^  -  P)^  +  2(E^  -  P)  (G  +  Pr^) 


-  ieA(l  +  F  -  r^A  )  (G  +  PrQ){ 


+  (G  +  Pr^) 


(1  +  F-r^A^)  (E    +  Pr^)  -  (E^  -  P)^ 


Solving  for  U(s)  using  Cramer's  rule  gives 


(B-4) 


U(s)     =    ^ 


%  ,  isd  -. 

-■:— (e      -1) 

IS  ^  ' 


-i(E^-P)s 


(G+F-rA^)s^-ieAs 

+Eo+Pro 


i(G+Pr^)s 


or,  after  expanding  the  determinant, 


U(s)     =    f 


,  isd  ,.  r 


(D  -IA^)s^ 

^    1     o    ' 


-i(G+Pr^)s 


(D  -I  A^)s^+(G+Pr  ) 

^    1     o     '        ^  o' 


(B-5) 


(G  +  Pr^)  (D^  -  I^A  ) 

+  (E^-P)(D^  -Io^^|s" 

+  (E^  -  P)  (G  +  Pr^)^ 


] 


(B-6) 
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Defining  the  coefficients  in  A  as 


c       =     (G  +  Pr  ) 

o  o 


(1  +  F  -  r  X^)  (E    +  Pr  )    -    (E    -  P)^ 


c^     =     -e\(G  +  Pr^)  (1  +  F  -  r^X^) 


(E^-P)    +    2(G  +  PrQ) 
+    (G  +  Pr^)  (1  +  F  -  r^r )  (F  -  Pr^  -  rX^) 


c^     =     -€Xm    -  I  \  )c 

3  ^     1  O       ' 


c^      =     (G  +  F  -  vX^)  (D^  -  I^A^)c 


c 

then 

A 
where 


=    1    +    F   -    D^    -    rX^    +   IX^ 

11  o 


s^  D(s) 


4.3  2  . 

Tiis)     =cs     +1CS      +CS      +ies   +   c 

^    '  4  3  2  10 

Using  this  nomenclature  Equation  (B-6)  can  be  written 
isd 


U(s)     _     1-  e' 


(D^  -  l^X^)  (G  +  Pr^  +  E^  -  P)s^ 


+    (G  +  Pr  )  (E     -  P) 

o  1 


(B-7) 


(B-8) 


(B-9) 


(B-10) 


Now  W(s)  is  determined.    Using  Cramer's  rule  again,  the  transformed 
radial  deflection  is  given  by 


W(s)     =    ^ 


(1  +  F-r  A^)s^ 
^  1     ' 


i(E^-P)s 


(D    -I  A^)s^ 

^1      o      ' 


0 


■r^(l-e^^^ 

IS  ^  ' 


(D    -I  X^)s^ 
^    1      o     ' 


-i(G  +  Pr^)s 


(D    -I  A^)s^  +  (G  +  Pr  ) 

^    1      o     '         ^  o' 


(B-11) 
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Expanding  the  determinant  gives 


W(s)     =    l^(i.e^«d) 


(1  +  F  -  r  A^)  (D    -  I  \^)s^ 
1  1        o     ' 


+    (1  +  F  -  r  \^)  (G  +  Pr  )s^ 

1  o' 


P.-Io^Ts^ 


(B-12) 


or 


W(s)     _    il 


isd. 
-^     1 


(D^  -  Iq^^)  cs^  +  (1  +  F  -  r^X^)  (G  +  Pr^) 


^o  isD(s) 

Finally,  the  rotation  expression  is  determined  as 


(B-13) 


Hs)    =    4- 


(l+F-r^\^)s^ 

-i(E^-P)s 

0 

i(E^-P)s 

(G+F-r\^)s^-ieXs+E  +Pr 

^                         '                              0            0 

IS  ^              ' 

(D  -I  A^)s^ 

^     1      o      ' 

i(G+Pr^)s 

0 

(B-14) 

or 


Hs)     =    - 


/I       isd. 


(E     -  P)  (D    -  I  A  ) 

^    1  '  ^    1        o 


+    (G  +  Pr^)  (1  +  F  -  r^A  ) 


(B-15) 


Ms)    ^    _  a-e'^^ 

%  D(s) 


(E^  -  P)  (D^  -  I^A  ) 


+    (G  +  Pr  )  (1  +  F  -  r  A  ) 


(B-16) 


APPENDIX  C 
EVALUATION  OF  THE  DISCRIMINANT  OF  A  FOURTH  ORDER  POLYNOMIAL 

th 


The  discriminant  of  an  n     order  polynomial 


ax+ax        +...+a       x  +  a       =     0 
o  1  n-i  n 


(C-1) 


is  defined  as  (15) 


A     =     a' 


2n-2 
b 


S  S 

n-i  n 


n-i 


2n-2 


(C-2) 


The  preceding  determinant  is  called  Vandermonde's  determinant  and  the  elements 
of  it  are  defined  as 

Sq     =     n  (C-3) 


(-4)' 


2a 


3a 


ka. 


\-x 


a  0 

o 


a         a 
1  o 


0 
0 


(k-  1,2, ...,n) 


(C-4) 


For  a  fourth  order  equation  of  the  form 


4  2 

ax+ax+a      =     0 

O  2  4 


(C-5) 
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the  S,  's  become 


'^  =    (-t)^ 


—  \  a       =     0 


(C-6) 


i-i)^ 


2a  0 

2 


2a  a  a 

a?       -         -o 


(C-7) 


(-t) 


0 
2a, 

0 


a  0 

o 

0  a 

o 

a_  0 


=     0 


(C-8) 


or 


0 
2a, 

0 
4a 


o 
0 


[4a^  a    -  2a    a^ 
o     4  o     2 

J 


0 
0 


2a 


o^ 

4 

0 

0 

\ 

O 

4a 

a„ 

0 

—  (2a    a    -  a^) 

2  O       4  2 


(C-9) 


(C-10) 


(-t) 


0 

a 

0 

0 

0 

0 

2a 

2 

a 
o 

0 

0 

2a 

0 

a 

0 

0 

2 

o 

1 

0 

0 

% 

0 

0 

a 

0 

a 

0 

2 

o 

4 

4a 

4 

a 

2 

0 

a 
o 

4a 

0 

a 

0 

a 

4 

2 

o 

0 

0 

^2 

0 

0 

a. 

0 

a„ 

0 

(C-11) 


or 
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2a 


1 

4a 

a 

a 

3 

4 

2 

o 

O 

0 

0 

0 

(0)     =      0 


(C-12) 


S,      =     — 


0 

a 
o 

0 

0 

0 

0 

2a 

2 

0 

a 
o 

0 

0 

0 

1_ 

0 

a 

2 

0 

a 

0 

0 

0 

6 

^O 

4a 

4 

0 

a 

2 

0 

% 

0 

0 

a 

4 

0 

^2 

0 

a 

0 

0 

0 

a. 

0 

a 

0 

5 

a 


2a 


0  0 


2 

o 

2a 

0 

0 

a 
o 

0 

0 

4a 

4a 

4 

s 

0 

a 
o 

0 

1 

4 

a 
o 

0 

0 

0 

a 

2 

0 

^o 

0 

0 

a 

0 

a 

0 

— 

o 
0 


0 

0 

a 
o 

0 


2a 

a 

0 

2 

0 

J_ 

4a 

a 

a 

_  J_ 

3 

4 

2 

o 

3 

a 

a 

o 

0 

0 

a 

4 

a 

2 

2a    r 

3a    a 

3           2 

0      4 

a 

0 

^- 

— 

2a^(a^  -  a    a  )  -  a  (4a    a  ) 

2^    2  O       4'  O  2       4' 


(C-13) 
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The  discriminant  can  now  be  evaluated. 


A     =     a 


4 

0 

-2^ 
a 
o 

0 

0 

a 
o 

0 

2                    2 
-    2(2aoa,-a,) 

2^^ 
^o 

0 

a 
o 

2 

0 

0 

2 

2 

a 
o 

(2a  a  ■ 

^        O     4 

-<) 

0 

2a 

— ^(a^-3a  a  ) 
?.      2       04' 

Expanding  the  determinant  about  the  last  column  gives 


(C-14) 


A     =     a 


—  (2a  a    -  a^) 
2        o  4        2 

a 
o 


a 

-2^ 
a 
o 


a 
-2^ 


-—(2a  a  -a^) 
2  ^     o  4     2' 
a 
o 


2  2 

-=-(2a  a  -a^) 

2    ^        0     4        2 

a 
o 


2^     2 

— ^(a    -  3a  a  ) 
3^  2         04 

a 
o 


.2!i 
a 
o 


-2!^ 

a 
o 


a 

a 
o 


—  (2a  a  -a^) 

2  O     4        2 

^O 


(C-15) 
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Expanding  this  further  gives 


A     =     a< (2a  a    -a) 

o  J        4  ^      o    4         2' 

^      o 


2  -I 

8/0  2,       .  ^2 

(2a  a    -  a  )  -  4  — 

2^04        2  2 

^  a  a 

■-    o  o  -■ 


4a 

+  — ^  (a^  -  3a  a  ) 
4^2  04' 

a 
o 


-—(2a  a    -  a^)  -4  — 
2  ^     o   4        2'  2 

a  a 

•—      o  o 


(C-16) 


A     =     -4  (2a  a    -  a^  ")       [-16a  a    +  8a^  -  4a^] 

\04  2/  04  2  2' 


+4a^(a^  -  3a  a  )  (-16a  a    +  8a^  -  4a^) 
2^  2  04'^         o  4  2  2' 


=     16(a^  -  4a  a  )   -(2a  a    -  a^)     +  a^(a^  -  3a 
2  o   4  o   ^         '^  '^     '^ 


4        2' 


2'  2 


a  ) 

O     4'J 


=     16(a    -  4a 

2 


.  a  )    -4a 
o  4  |_      o 


2  2  44 

a    +4aaa    -a    +a    -3a 

4  0    2    4  2  2 


.  a  a 
o  2  4J 


2  222 

=     16(a    -  4a  a  )  (-4a  a    +  a  a  a  ) 
^2  04'^       04        024' 


2  2 

16(a    -  4a  a  )  (a  a  )  (a    -  4a  a  ) 
^2         o  4'  ^  o  4'  ^  2         o  4 


Finally 


—  22 

A    =     16a  a  (a    -  4a  a  ) 
o  4^  2         04 


(C-17) 


APPENDIX  D 
DETERMINATION  OF  THE  CRITICAL  VELOCITY  EQUATIONS 

Equation  (4-3)  is  the  first  expression  satisfying  the  conditions  for  which  re- 
peated roots  occur .    It  is  repeated  here  for  reference . 


c     -  4c    c 

2  O       4 


0 


(D-1) 


With  the  definition  of  terms  as  given  by  Equation  {D-2), 


e       =     E    +  Pr 
o  o  o 


^     =     E^-P 


e       =     E     -  P  +  2(G  +  Pr  ) 

2  1  ^  O' 


e^     =     G  +  Pr 

3  o 


e       =     I    -  r 

4  o  1 


e       -     e  f    -  e  e 

5  O   1  12 


e       =     e  f    -  e 

6  o  o        1 


=     1  +  F 


=     1  +  F  -  D 
1  1 


=     F  -  Pr 
2  o 


=     G  +  F 


=     D  r     +  I  f 
4  1    1        o  o 


r  f    +rf 

12  O 


-     e  D    -If 

6  4     10   1 


7 

— 

1^1  -  V 

e 

= 

f    -rx^ 

O           1 

9 

- 

D  e     -  I  f 

14            0    1 

(D-2) 


111 
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The  coefficients  in  Equation  (3-10)  can  be  written  in  the  following  form. 

c       -     f    +  e  X^ 

1  4 

2 

c       =     e  (f  e    -  e  ) 

0  3    e  o        i' 

c       =     -eXe  f 

1  3  8 

c       =     f(ec-ee)+ef(f    -  rA^) 

2  Y      O  12  3  e     2 

c       =     -e\f  c 

3  7 

c       -     (f    -  rA^)f  c 

4  3  7 

c       =     f  (e    +  e  ) 

5  7       3  1 

c,     =     e  e 

6  13 

c       =     f  (f    +  e  A^) 

7  7^    1  4        ' 

c       =     e  f 

e  3  8 

c       =     e  f    +  c  {D-3) 

9  17  8 

The  coefficients  in  Equation  (D-3)  are  substituted  into  Equation  (D-1)  and 

2 

terms  are  collected  which  multiply  various  powers  of  A   .    Proceeding  to  this 
end,  the  first  coefficient  is  c    . 

2 

c^     -     (D    -  I  A^)  [e  (f    +  e  A^)  -  e  e  ]   +  e  (f^  -  r  A^)  (f    -  rA^) 

2  ^1  O'O^l  4  12  30  12 

=     Def    +DeeA     -Dee 

lOl  104  112 

-  I  e  f  A^ 
o   O  1 

2  4 

+  IeeA  -leeA 

o    1    2  o    O    4 

+  e  f  f    -  e  (r  f    +  rf  )A^      +  e  rr  A^ 

302  312  O  31 
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c       =     D(ef-ee)+eff 

2  lOl  12  3O2 


+  (Dee    +  I  (e  e     -  e  f  )  -  e  (r  f    +  rf  )lx' 

104  O^     1    2  01'  312  O" 


+  (e  rr     -  I  e  e  )A 
^31       004' 


c       =     De    +eff    +(Dee    -le    -ef)X    +ferr    -lepu"* 
2  15        302^104       o  5        35'^        ^31       00    4^ 


c^     =     e    +eX    +e\"* 

2  y        e  9 


where 


e       =De+eff 

7  15  3  O  2 


8      =     Dee    -le    -ef 

e  104       05        35 


e       =     err-Iee  m-4^ 

9  3     1       o  o   4  \^  ■*> 

Now  this  coefficient  is  squared 

cl     =     S  +  2e^e^A^  +  (ej  +  2e^eg)x'' +  2e^eg\'' +  e^A^^  (D-5) 


c       =     e[e(f-r\)-e] 
o  3^  o^  o        1     '        1' 


/     r         2,  2 

=     e^(e  f    -e)-eerA 

3      O  O  1  3    0    1 


=    ^<\-v/)  (D-6) 


c       =     (f    -  rx^)  (D    -  I  A^)  (f    +  e  A^) 

4  ^3  '^1  0^1  4 


=     (f    -  rA^)  [D  f    +  (D  e    -  I  f  )A^  -  I  e  a''] 

3  1114  01'  04 


=     f  D  f    +  [f  (D  e     -  I  f  )  -  rD  f  ]A^ 

311  '3^      14  01'  1    1' 


+  [r(I  f    -  D  e  )  -  f  I  e  ]  A^  +  rl  e  A^  (D-7) 

Oil    4304'  04  ^  ' 
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Let 


D  e    -  I  f       =     f 

1    4  O   1  9 


then 


4  6 

; 

4 


c       =     f  D  f    +  (f  f    -  rD  f  )X    -  (rf    +  f  I  e  )A     +  rl  e  A 

4  3      1    1  '    3   9  1    l'  ^9  3  0     4  O     4 

and 


4c  c       =     4e    {(e    -  e  r  X  )c  } 

O    4  3  ■■      6  O    1       '    4-' 

=     4e^  {(e^  f  D  f  )  +  [e  (f  f    -  rD  f  )  -  e  r  f  D  f  ]  x^ 

3  6  3     11'         ^    6^  3  9  1   l'  O    1   3     1   l' 

-  [e^(rf    +  f^I  e  )  +  e  r  (f  f    -  rD  f  )]  a^ 

6^      9  3  0    4'  O    1^  3  9  1   l'' 

+  [e^rl  e    +  e^r  (rf    +  f  I  e  )]A^  -  e  r  rl  e  A®  } 

6      0    4  O    1^      9  3  0    4'  O    1      O    4        J 


Finally 


(D-8) 


2  4  6  e 

4c  c       =e       +eA+eA+eA+eA  (D-9) 

0    4  lO  11  12  13  14  ^  ' 


where 

e  =     4e  e  f  D  f 

10  3    6  3     11 

e  =     4e  [e  (f  f    -  rD  f  )  -  e  r  f  D  f  1 

11  3'    6'   3  9  11'           O    1   3     1   1' 

e  =     -4e   [e  (rf    +  f  I  e  )  +  e  r  (f  f    -  rD  f  )1 

12  3      6'      9          3  O    4             O    1^  3  9  1   1   ' 

e  =     4e   [e  rl  e    +  e  r  (rf    +  f  I  e  )] 

13  3^    6      O    4  O     1^       9  304'' 

e  =     -4e  e  r  rl  e                                                                                        m-10) 

14  3  o  1     o  4  \^  ^"J) 

Equation  (D-1)  now  becomes 

C  A® +  C  A^  +C  A^  +C  X^  +C       =     0  (D-11) 

6  6                   4                  2                  0                                                                                                                   ^               ' 
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where 


C       =     e^  -  e 

e  9        14 


C       =     2e  e     -  e 

6  8    9  13 


C       =     e^  +  2e  e     -  e 

4  8  7     9  12 


C       =     2e  e     -  e 

2  7  e        11 


C       =     e    -  e  (D-12) 

O  7  lO  ^  ' 

If  the  axial  and  rotatory  inertia  are  neglected 


r 

1 

0 

^o 

= 

0 

^ 

= 

0 

K 

= 

0 

f 

4 

= 

0 

^5 

= 

^^o 

The  coefficients  in  Equation  (D-11)  reduce  to 
C^     =     C       =0 

6  8 

C    -  C       =     e   [e    r^  f^]       =     e^  r^  f^ 

4  4  3       3  O  3  O 


^2 


(D-13) 


C       =     2  {D  (-e    e    f  r)  +  e   [2e    r  D    f    -  e    f    f    r]} 

2  1350  36  11  SOs'-" 

C    -*  C       =     2e   [2e    r  D    f    -  f    r(D    e    +  e    f    f  )) 

2  2  36  11015  30     2 

C       =C       =(De+eff)^-4efDfe  (D-14) 

O  O  1530     2  33116 


APPENDIX  E 

COMPUTER  PROGRAM  FOR  SOLUTION  OF  LOAD  VELOCITIES 

WHICH  CAUSE  REPEATED  ROOTS  IN  THE  UNDAMPED 

CHARACTERISTIC  EQUATION 


The  program  shown  here  is  used  for  calculation  of  the  location  of  the  re- 
peated roots  of  the  undamped  characteristic  equation  using  the  polynomials 
defined  in  Chapter  IV.    The  program  is  written  in  Fortran  for  use  on  General 
Electric  Company  Mark  II  time  sharing  system.    The  only  requirement  for  run- 
ning the  program  is  that  a  data  file  which  is  named  "EMFILE"  must  be  estab- 
lished which  includes  the  data  for  the  cases  to  be  solved. 

EMFILE  is  set  up  such  that  there  are  seven  numbers  in  each  line  of  the 
file,  the  first  being  the  file  line  number;  the  next  six  are  the  dimensionless 
parameters  h/R,  Eg  /Ex^,  G^^,  /-^Xq'  ^u  /^Xq'  -^'  ^^^  ^  ^^  *^^^  order. 

2 

Starting  the  program  will  cause  the  eight  roots  for  A    to  be  calculated  and 
printed.    Also,  the  program  automatically  calculates  the  set  of  roots  in  which 
the  axial  and  rotatory  inertia  are  neglected.    The  following  is  a  listing  of  the 
program . 
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v/C'<ir 

100    DIMENSION    A(5),n;'^(4),C''?C4) 

110     b    KEAOC"EMFILE"*J?)DUMy/»H>ErOf:X*G^ZP1F':X,R:i  ILP^F 

120    2    FO'-^MATCV/) 

I 

130  ^=4 

140  IF(P-0)7,8,8 

150  7    RO  =  H/2;G0T-3     9 

160  H    ►^0  =  -H/2 

170  9    01  =  CH)**2/12 

IHO  XKX=    3. 1 4159/(12.  )t. 5 

190  n=(XKX)T2*GXZ0EX 

200  PRI>JT*"  H  ETHETA/EX       GXZ/EX  E^JU/EX  P 

210  &  F  G" 

220 

230  PRIMT    lO,M,ET0EX,GXZ0EX,Filj,p,F,G 

240  10    F0r<MAT     (7F10.5) 

250  PRINT*"*********** ********♦♦♦*♦* ♦*♦♦******♦*♦******* 

255  &*****************'• 

The  above  steps  read  the  data  from  EMFILE  and  print  the  data  used  as  a 
heading  on  the  output. 
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260  P^IMT* 

270  EOIJ=ET0EX*C1+HT2/12) 

280  R=H 

290  XI0=Ht3/12 

300  R1=H 

310  EO=EOU+P*RO 

320  E1=E1U-P 

330  E2  =  EUJ-P  +  2*(G+P*R0) 

340  E3=G+P*R0 

3b0  15    E4=XI0-R1 

360  Fl=l+F-01 

370  F2=F-P*R0 

380  F0=1+F 

390  F3=G+F 

400  E5=E0*F1-E1 *E2 

410  E6=E0*F0-E1 *E1 

420  F4=D1*R1+XI0*F0 

430  F5=R1*F2+R*F0 

440  F6=E4*D1-XI0*F1 

450  F9=D1*E4-XI0*F1 

460  E7=D1*E5+E3*F0*F2 

470  E-8=D1*E0*E4-XI0*E5-E3*F5 

480  E9=E3*R*R1-XI0*E0*E4 

490  E10=4*E3*E6*F3*01 *F1 

500  El 1=4*E3*(E6*(F3*F9-R*D1*F1 )-E0*Rl*F3*Dl *F1 ) 

510  E12  =  -4*E3*(E6*(R*F9  +  F3*XI0*E4)+E0*R1 * ( F3*F9-R*D  1  *F  1  ) ) 

520  E13=4*E3*(E6*R*XI0*E4+E0*R1*(R*F9+F3*XI0*E4)) 
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530    H:1  4  =  -4*E3*F:0*R*R1*XI0*E4 

bAQ    C0=E7+E7-Ein 

550    C2=2*E7*E8-E1 1 

560    C4=E8*E8+2*E7*E9-E12 

570    C6=??*E8*E9-E13;    C8=E9*E9-El/i 

580    C0P  =  F1  *F3+n>l 

590    C2P=E4*F3*D1-F1 ♦CR*D1 +F3*X 10 ) 

600    C4P=F1*R*XI0-E4*(R*D1 +F3*XI0) 

610    C6P=E4*XI0*R 

630    PKIMT    82 

630    82    FOR.^ArC/"  CO  '  C2  C4 

640    4  C6  C8"     ) 

650    90    PRINT    20.C0#C2,C4, C6,C8 

660    20    F3RIMATC5E1  4.  6//) 

The  above  steps  calculate  the  coefficients  for  the  polynomials  and  print 
those  for  the  fourth  order  polynomial , 
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'S70  Ad  )  =  nojA(y )  =  r:  j;A(3)=n4;A(<Q)  =  n6;A(S)=i; '^ 

690    P'<T'\)T, 

700    PRINT,"  REAL    COMH^^J^lM  I 

705    5IMAGIMARY    COMPOMEiNJT" 

710    00    2  5    I=1,M 

7^0    ?.5    P^INT    50,  I,  ^^(  I),  C^^  CI) 

73:)    50    FOXMAI     (5H    -M,-)  T,  !;:>  >  1  0  X,  1  P2E20  .  7  ) 

740    PS'InJT    60 

750    60    FORMA  I  C/"     IHE    ^^EC^'MST.     COEFF.     OF    THE    POLYM.     A'^Ef) 

760    Pf^INl, 

770    /i  =  ,M+l 

780    L=M+1 

790     IF((-1 )*+M)l 10, 190, 100 

800     too    Pi^l,\]r     102 

810     102    FOi^MAK/"  CO  C2  C4 

B20    &  C6  C8") 

The  above  steps  solve  for  the  roots  of  the  fourth  order  polynomial  and  print 
the  reconstructed  coefficients  as  a  check  on  the  validity  of  the  roots. 
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830     IF     (-^-:>)67,  73,  73 

HAO    67    DO     70     I=L*5 

RSO    70    A(I)=0 

860    73    P-^INT    20,A(  1  ),  A(2),A  C3),AC4),A(',>) 

870    P.r^IIMT, 

880    i^^IXT     104 

JtgO     li)/i    Ffl^MAfC  COP^TME  C2P'UMF:  r./^prs'iv,:-: 

895    *C6P'>ITMF") 

900    PJUMT    66,C0P,C2P,C4P,C6P 

910     66    H"j;>;  ■iAr(4Fl  4.  6// ) 

9;£0    A(  1  )=C0F;  A(2)=C2PJA(3  )  =  C4P;A{4)=C6P 

930    7S    N=i\)-1 

940     IF(2-N)80, 130* 1  40 

9b0     140     IF(-,M)80,  1  50,  1  bO 

960     110    PF^INJT     112 

970     112    F0RMAr(/"  COPRIME  C2PRIME  C4PRIMF 

975    &C6PRIME") 

980    IF(M-5)72,83,83 

990    72    0081     I=L,5 

1000    81    ACI)=0 

1010    83    PRINT    20,A(1  ),A(2)*A<3),A(4),A(5) 

1020    PRINT, 

1030    G0T0    75 

1040    130    VCR  =  E6/(p:0*R1  ) 

1050    PRINT    120,VCR 

1060    PRI^JT    122 

Lines  830  to  1060  direct  the  solution  and  printing  of  the  other  four  roots 
for  ^^ . 


122 

1070  IP.P.    FGRMATC    ,\)EGLECT.     R0TAT.     Af^D    AX.      I-nJE^TIA,      THE 

1080  &    CTF.FF.     ".    f^00TS    ARF:"/) 

1090  PRINT    99 

1100  99    F?)KMAT("  *****  +  ***>4:  +  ****:>  +  .i-*  +  -+-^ 

1105  5+***+*+") 

1110  IFf-XI0)30,40,40 

1120  30    CONTIiNJUE 

1130  XIO=0 

1140  R1=0 

11  SO  GOTO     15      . 

1160  1     )^RI.NT> 

1170  150    PRINT* 

1180  190    FOR.^ArC    THE    FINAL    CRITICAL    VELOCITY    IS    "F15.8//) 

1190  PRI^JT    152 

1200  152    F'5RMAT("    THE    FINAL    GRIT.     V/EL.    APPR.     INFINITY"//) 

1210  PRI^]T    160 

1220  160    F0R.V1ATC" 

1230  & ••///) 

1240 

1250  40    G0T0     5 

1260  300    ST0P 

1270  FNO  ' 

This  section  of  the  program  makes  the  simplification  of  neglecting  axial 
and  rotatory  inertia  and  directs  the  solution  for  the  corresponding  roots . 


APPENDIX  F 
PARTIAL  FRACTION  EXPANSION  OF  A  FOURTH  ORDER  POLYNOMIAL 

The  characteristic  equation  D(s)  can  be  factored  as 

D(S)     =     c^(s  -  s^)  (s  -  s^)  (s  -  S3)  (s  -  s^)  (F-1) 

so  that 


(F-2) 


2  2 

c's+c'  cs+c,  ... 

5 6       _       5 6 _       A(S) 

^(s)  "     (s  -  s^)  (s  -  Sg)  (s  -  S3)  (s  -  s^)     "     B(s) 

Let 

B(s)  s  -  s^         ^  "  ^2         s  ~  S3         s  "  ^4 

The  coefficients  j3j^  can  be  found  (16)  from  the  expression 

A(Sk) 
/3.       =     (F-4) 

B'(Si^) 

which  means  that  A  over  the  derivative  of  B  is  to  be  evaluated  at  s  =  Sj^.     This 
derivative  is 

d    „,  >  ,  .   d 


B 


(s)     =    -^  B(s)     =     (s  -  s^)  ^  [  (s  -  s^)  (s  -  Sg)  (s  -  s^)  ] 


+    (s  -  s^)  (s  -  S3)  (s  -  s^)  (F-5) 


Furthermore 


^[(S-S^)(S-S3)(S-SJ]        =       (s-S^)^[(S-S3)(S-S^)] 

+    (S  -  S3)  (S  -  s^) 

and 

^[(S  -  S3)  (s  -  s^)]      =     (s  -  S3)    +    (s  -  s^) 
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Therefore 

d    ,, 
^  [  (s  -  s^)  (s  -  S3)  (s  -  s  J]      =     (s  -  s^)  [  (s  -  S3)    +    (s  -  s  J  ] 

+    (s  -  S3)  (s  -  s  J 
and 

B'  (s)     =     (s  -  sj   ^(s  -  s^)  [  (s  -  S3)    +    (s  -  s  J]  +    (s  -  S3)  (s  -  s^) 
+    (S  -  S^)  (S  -  S3)   (s  -  s^) 
Finally  the  coefficients  can  be  condensed  to  the  expression 

^k     -     ~  -     -^ .  P    7^   k  (F-6) 

P=l 


APPENDIX  G 

CHECK  TO  SEE  THAT  SOLUTIONS  SATISFY  THE  GOVERNING 
DIFFERENTIAL  EQUATIONS 


The  governing  differential  equations  have  the  form 


f    dfu^e    ^  +  f   ^    -   0 


(G-1) 


dU     ,f  2   dW     ^,  dW  ^      „, 

e    -jT  -  (i    -  rX  ) ex  -rr  +  e  W  -  e 

1  d0      ^  3  '  ^^2  d0        o 


3^    =   -qjH(-0)-H[-(04-d)] 


-,    d^U  dW  ^  ^    d^ 

i_ -  e„  -r-r  +  f^  — ^  - 


■"dc/,^ 


d<p 


^  dcp^ 


e^ip  =  0 


(G-2) 
(G-3) 


First  find  the  second  derivatives, 
k 


K  (<^)       =     i-\  +  ''\f  Fj^(0)     =     -s^  F^{4>) 


K  i<t>)    =    *^— ^  F,  (0) 

w  k  ^^'  CO!,         k^^' 

q    k 


u^k  <'^)   =  r  \^^^ 

9 


(G-4) 


Substituting  into  the  first 


8  C 


E-sfF,  ((^)     +E  — F,    ^+e    E— ^F,  -f    Z)  sf  F, 
,         k     k^^'        1    c„     k  J        1  ,    c„a,      k      7  V    k     k 
k  -ik9-^  kak  k 


? 

=     0 


4 

E 

k=l 


c  f  c 

8  c       k       c 

9  9 


TT  (^-  ^l  +  °  J    -    t  ^f 


C_    '    7      k  8  7      k 


Fk(0) 


? 

=     0 


The  bracketed  term  must  be  zero  if  Equation  (G-3)  is  satisfied  by  the  solution. 
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(c  e  c  \ 

f   -^    -    ^^    -    f   1 
« %       ^9       7 


f  c    -  e  c       o 

.        6     6  IB        ?        - 

s,     +  =     0 


(G-5) 


Look  at  the  second  term  to  see  if  it  is  zero  because  both  terms  in  Equa- 
tion (G-5)  must  be  zero. 

fc     -ec       =     fc,    -eef      =     f(c,-e  c„) 

86  IB  86  138  8^6  13' 

=     f   (e  e^  -  e  e  )     =     0 

8^13  is' 

Now  the  first  term  is  inspected. 


f„c    -  e  c    -  f  c 

8     5  17  7     9 


f   [f   (e    +  e  )  -  e    (f    +  e  A  )  -  (e  f    +  e,f  )] 

7^8^3  l'  1^1  4        '         ^17  3   8" 


=    f   [e   (f    -f    -  f    -  e  \  )] 

7        18         7  1  4       " 


=     ef[f    -rX^-D    +IX^-1-F+D 

1  7       O  1  10  1 


-  (I    -  r  )\  ]      =0 
^  o       i'     ' 

Thus  the  solution  satisfies  Equation  (G-1).    Now  Equation  (G-3)  is  checked. 


4 

E 

k=l 


c^              f^c,  (c    s,   +  c  ) 

o52,76  ^7k           8'            r        2 

1^  —  s,    +  +  e„ -   f    s,     -    e 

7                                         Cg                                                                                           k                                                                                                                                                                                                                                Cg  3                                                                                                                                                                                                                                  Cg                                                                                                                                                                                                                                                                                                                                                                                                                                  7                                                                                         k                                                                                                                                                                                                                                  3 


(G-6) 


or 


(c  e  c  \    „        /f 

7  Cg  Cg  7/     k  V     ^ 


^-^6  ^3^8                   \        ?      , 

-  +  -    e    1    =    0 

c  c                3, 

9  9 


,)    = 


Look  at  the  second  term. 

fc     +e(c-c)     =     e(fe+c-c)     =     e(ef-ef)     =     0 

76  38  9  3^71  8  9'  317  17 

Now  look  at  the  first  term 


f   [c     +   e   (f    +  e  \  ) 

7         5  3      1  4        ' 


-ef     -ef]     =f[fe     +e(f+e\-f+f)-fej     =0 

17  3   8  7       7    1  3^  1  4  8  7  7    1 
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So  Equation  (G-3)  is  satisfied  by  the  solution.     Now  Equation  (G-2)  will  be 
checked. 


c       —    k^  k       ^  (     CgSj^  ^k  "*■  ^k  ) 


^  i  E  s  F     +  E  ( - 

'9        k       '^    '^  k    \     "g'k 


(f    -  r\  ) 

^  3  ' 


is^(c    s.   +  c) 


_  Jk^  7  '^k  '  "a' 
k  9 


^-       "     'c"        °      ^k 


-     €Xl] 


{-•■^'.)-'. 


i(c  s,    +  c  ) 
^  7  k       s' 

F,     +    A, 

s,  c  k  k 

k  9 


e    E(i\F^     =     -    {H(-0)    -    H[-(0+d)]} 


(G-7) 


where 


sgn  (b^  c^aj^ 


2  s, 


{H[sgn(bj^)0]    -    H[-sgn(b^)0] 


-    H[sgn(bj^)  (0  +  d)]    +   H[-sgn(b^)  (0  +  d)]} 


/3 


Aj^     =     sgn(bj^   2i^  {-H[sgn(bj^)0]     +   H  [ -sgn  (bj^)c^  ] 


+   H  [  sgn  (b^)  (0  +  d)  ]    -    H[-sgn(bj^)  ((/)  +  d)]} 


Collecting  terms  gives 


k=l 


I — le.  c   s 


,  ^   ^,  e  c,        f    -  r\  ) 

15k  .16  3  ' 

-    1 + 


c  s, 

9    k  9 


is,   (c    s,   +  c  ) 
c  k  ^  7    k        e 


^  le              ., 

+  —  (c    s,     +  c  )    +  (c    s,     +  c  ) 

c    ^  7    k  e'         s,  c    ^   7  k            ' 

9  k  9 


ie    s, 

3   k 


^k   ^   W   ^   %\ 


{H(-0)    -    H[-((^+d)]} 


(G-8) 
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Looking  first  at  the  term  multiplying  F.  ,  which  is 

*-     9 


+   e. 


eXc  eXc 


s,    +  — 

c         k  c 


(e  c    -  e  c  J 

8      .      .    ^    O    8           16' 
—    +     1  


Multiply  and  divide  this  by  (is, )  and  get 

— I—  {  -  c  (f    -  rX^)  sf    +    ie\c    sf   +    sf  [  e  c    -  c    (f    -  rA^) 

CIS,       '^  73  k  7k  k'l5  83 


e^c    +  e    c 

0    7  3      9 


+   ieXc    s,     -    (e  c    -  e  c  ) } 

8    k         ^  o  a        1    6'  •■ 


Look  at  the  coefficients  on  powers  of  s. 


-  e  c     +   e  c, 

O    8  16 


e   (e  f  )    +  e    (e  ej 

O  ^    3  8'  1^13' 


e(fe    -e)     =    -c 

3^83         r  o 


e\c„    =    e  f  eX   =     -c 

8  3   8  1 


e  c     -   c    (f    -  rA.  ) 

15  8^3  ' 


-   e  c     +   e  c 

O    7  3    9 


=    fe(e+e)-ef(f-rX)-    e  f  c 

71^3  l'  38^3  '  OY 


+   e    (e  f    +  e  f  ) 

3       17  3  8 


-f   (ec-ee)    -   ef   (-e    +f    -  rX  ) 

70  13  383  3 


-  f^  (-  e  6    -  e  ) 

7  ^        13  l' 


f   (e  c  -  2e  6    -  6^)    -   e  f   (-G  -  Pr  +  G  +  F  -  r\^ ) 

7^0  13  1  38  O 


f   [e  c  -  e    (6    +2e  )]    -   e  f   (f    -  r\  ) 

70  1^1  3'  38^2  ' 


-f   (ec-ee)    -ef(f-r\)     =     -c 

70  12  3  8      2  2 
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Similarly 


eXc  =     -c 

7  3 


c    (f    -  rA^)     =     -  c 


SO  the  term  multiplying  F,   is  the  characteristic  equation  which  is  evaluated  at 

s  =  s,  .     Of  course  this  equation  goes  to  zero  for  all  roots  s,  . 

The  remaining  identity  which  has  to  be  proven  is 

4 

E  (e^Bj^  +  c^Aj^)     =     -{H(-0)-H[-(0+d)]}  (G-9) 

k=l 

Noting  the  identity 

sgn  (b  ) 

2"-^  {  H  (  sgn  (bj^)  (^  ]    -   H  [ -sgn  (bk)  (^  ]    -    H  [  sgn  (bk)  (0  +  d)  ] 

+    H[-sgn(bk)  ((^  +  d)]}     =     -   {H(-0)  -  H[-(0  +d))}        (G-10) 
which  holds  whether  b,   is  less  than  or  greater  than  zero,  the  remaining  state- 


ment which  is  to  be  proved  is 

y^   /_i   6    k  o*^  k 

k=l  \ 


s,  s, 

k  k 


=     1 


or 

4     a 
Z-«   —     ec      -e(cs,   +c)        =     1 

Substituting  for  c.  (i  =  6,  7,  8)  gives 

4     a, 
y;   -^  [  e  f   (f    +  e  A  )  sf  +  c    ]      =     -  c  (G-11) 

.*-'       S,07^14'k0  4  ^ 

k=l      k 
This  is  an  identity  and  is  programmed  to  be  checked  with  each  case  calculated 
using  the  computer  program. 


APPENDIX  H 
COMPUTER  PROGRAM  FOR  DEFLECTION  AND  STRESS  CALCULATIONS 

The  computer  program  DEFSTR  calculates  the  deflections  and  stresses  in 
the  cylinder.    There  are  two  required  functions  to  perform  in  order  to  operate 
the  program  as  set  up  for  operation  on  the  General  Electric  Company  Mark  II 
Fortran  time  sharing  system.    The  first  is  the  establishment  of  a  data  file 
(called  EMDATA  in  the  program) .    This  file  has  ten  elements  on  each  line  of 
the  file.    The  first  element  is  the  line  number  and  following  that  are  the  dimen- 
sionlessdatah/R,  Eg^/E^^,  G^z^E^^,  E^^/E^,  P,  F,  e ,  X^,  and  d. 

When  the  data  file  has  been  established  and  the  program  DEFSTR  called 
and  asked  to  run,  a  request  for  data  is  received  which  must  be  in  the  form: 
KEYl,  KEY2.    KEYl  and  KEY2  are  integer  values  which  allow  the  selection 
of  the  type  of  analysis  to  be  done  and  the  type  of  output  desired.    KEYl  deals 
with  the  inclusion  or  neglection  of  inertia  terms ;  KEY2  allows  the  amount  of 
analysis  and  output  to  be  regulated  as  shown  in  Table  H-1. 
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Table  H-1.    Options  Available  for  Program  DEFSTR 


Mode  of  Operation  of  the  Computer  Program 

KEYl 

Neither  Axial  Inertia  nor  Rotatory  Inertia  is  Neglected 

Both  Axial  Inertia  and  Rotatory  Inertia  Are  Neglected 

The  Axial  Inertia  Is  Neglected 

The  Rotatory  Inertia  Is  Neglected 

Neither  and  Then  Both  Inertias  Are  Neglected 

Neither  and  Then  Both  Inertias  Are  Neglected 

Neither  and  Then  the  Rotatory  Inertia  Is  Neglected 

Neither ,  Then  Axial ,  Then  Rotatory 

Neither,  Then  Axial,  Then  Rotatory,  Then  Both 

1 
2 
3 
4 
5 
6 
7 
8 
9 

KEY2 

1 
2 
3 

4 
5 
6 

7 

8 

The  Deflections  Are  Printed 

Deflections  and  the  Individual  Stresses  Are  Printed 

Deflections ,  Individual  and  Surface  Stresses  Printed 

Deflections  and  Surface  Stresses  Printed 

Only  Surface  and  Individual  Stresses  Printed 

Only  Individual  Stresses  Printed 

Only  Surface  Stresses  Printed 

Only  the  Roots  of  the  Characteristic  Equation  Printed 

A  listing  of  the  program  is  included  here. 


100    FilNlOTDJ    ■■1iJLr<^X,-^) 
110    -iiAL     -IMLT 
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120  IF(AX.EQ.O)G0T0     7000  J  IF  CB.  EQ.  n  )G'^T'^     7000 

130  AL=AL0G(ABS(AX) ) J gL=AL0G (A8S (B ) ) 

140  ALIM=-80.0 

150  IF(AL+gL.LT.ALIM)G0T0    7000 

160  Pi<0D  =  AX*9;G0  ra     7010 

170  7000    Pf^TDsO 

180  7010    MULT=PR3D 

190   i^etiir^^;e>jd 

yOO  FtJNCTISN    DI\/(AX,B,C*D) 

PIO  REAL    MULT 

P.P.n  Yl=v)iJLT(AX>C)  J  Y2  =  .v|ULTCB*D) 

230  CL=Y1+Y2 

240  Y3  =  f>^ULTCC#C);  Y4  =  Ml)LT  CD#  D) 

250  DL=r'3  +  Y4 

260  IF(CL.E-3.O)G3T0    7015 

270  AL=ALOGCABSCCL)) 

280  701 5  AL=0 

290  8L=AL3G(ABS(DL>) 

300  ALI'>1=-80.0 

310  IF(AL-3L.LT.ALIiv|)G0T0  7020 

320  QU0T=CL/DLiG0T0  7030 

330  7020  aU0T=O 

340  7030  DIV=QU0T 

350  RETURiMJEND 

The  steps  in  lines  100-350  define  two  functions  for  the  multiplication  and 
division  of  complex  numbers. 
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360  '^EAL     'lULF 

370  DIMEiMSI0N    A  (  5  )  ,  A  I  (  5  )  ,  R>?  (  4  )  ,CI^  C  4  )  ,  D^R  (3  )  ,  DC^  (3  )  ,  RA  (  4  )  , 

380  &CBC4),0  C121  ),EP1(?(4),E:P1  I  (  4  )  *  E:P2R  (  4  )  ,  EPf?  I  (4),SGi\(4), 

390  *FP2r^C4)*FP2I(4),FPRC4),FPI  (  4  )  ,  FLR  (  4  )  ,  FL  I  <4),FR(4),FI(4), 

400  <?.WFP1  I  C4),  WFP2R(4),  WFP2I  (  4  )  ,  WFPR  (  4  )  >  WFP  I  (4),  WFR(4^> 

410  ADFI  (  4),  DWFRC  4  ),r)WF  I  (4),  0?:P1RC4),I)RP1  I  (  4  )  ,  0FP2R  f  4  )  , 

420  &DUFR(4),0UFI  <  4  )  ,  UFPl  R  (  4  )  *  UFP  1  I  (  4  )  ,  IJFP2R  (  4  )  ,  UFP2  I  (  4  )  , 

430  ^    IIFPI  (4),IJFLR(4),UFLI  (  4  )  ,P'53R  (  1  2  1  )  ,  PI  ■)  T  CI  2  1  )  ,  WOOR  (  1  2  1  ), 

440  &UgOR<121  ),U30I  (121  )*SXVJ(121  ),ST'^(i;^l  )*SXg(121  )*ST9C121  ), 

4b0  ^>S  X:^  (  1  2  1  )  ,  S  r0  (  1  2  1  )  ,  S  X I  (  1  2  1  )  ,  S  T  I  (  1  2  1  ) 

460  ^.*'JFR(4)>UFI  (4),RR!  (4),CR1  (  4  )  ,  ■RR2  (  4  )  ,CR2  (  4  ) 

46  b  '<    CA(4),R=^(4),FP1R(4),FP1  I  (  4)  ,  WFPl  R  (4  )  >  WF  I  (  4  )  ,  DFR  (  4  )  , 

467  *    DFP2I (4)>UFPR(4), WIQI  (121  )>SX^(121  ) 

470  I\'PiJT,KFYl  *KF:Y2 

480  5    REAOC'EMDATA"*  1  0  )OUM\/,  H*  ET1  EX,  GXZ0EX,  E  1  ll,P,F,EP,V*0 

490  10    F0R'-1AT(V/) 

500  PRI>JT     15 

510  15    F0RMAT("  H       ETHETA/F.X       GXZ/FX       E^JU/^X  P 

520  ft  F  EPSIL0:M       LAM'3l.)AtP       D" ) 

530  PRIMT    2O>H,RTOFX,GXZ0FX,F1U*P>F,EP,  V/,D 

540  20    F0RMAr(F6. 4, 5F9. 6,F8.6,F8.3,F5.2) 

550  PRI'MT    30 

560  30    FORMAT ("*+***♦****♦***************♦****♦***+♦*** 

57 0  fiM*^.* ********************** ****"//  ) 

The  above  steps  reserve  the  required  locations  for  the  program  variables 
and  request  for  the  input  KEYl,  KEY2.    Then  the  first  set  of  data  is  printed 
and  calculations  begin. 
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b80  N=4 

590  L3=0 

600  IF(P-0)7,B.8 

AlO  7  f^0  =  H/2iG3T'3  9 

620  8  R0  =  -H/2  ^ 

630  9  D1=(H**2)/12 

640  XKX=  3. 1 41 59/(12.  )t. 5 

650  G=(XKX) t2*GXZ0EX 

660  E:0lJ  =  ETf^EX*Cl+HT2/12) 

670  R=H 

680  L=0 

690  G0r'^(4O,46#44*42»4O#4O*4O,4O*4O)»KEYl 

700  40    XI0=CH**3)/12 

710  R1=H 

720  0010    48 

730  42    PRINT    41 

740  41     F0RMAT<"    NEGLECTING    R0TAT0RY    INERTIA    GIVES;"//) 

750  XI0=0 

760  R1=H 

770  G0T0    48 

ISO  AA    PRINT    43 

790  43    F3RMAT("    NEGLECTING    AXIAL     INERTIA    GIVES:"//) 

330  XI0=CH**3)/12 

810  R1=0 

820  G0T0    48 

830  46    PRINT    47    ' 

840  47    F0RMAT("    NEGLECTING    AXIAL    AND    R0T.     INERTIA    GIVES:"//) 
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HbO    XI0=O 
860    ^\=n 
H70    G<?iT3     AB 
8P?0    48    L=L+1 

The  logic  in  lines  580-880  establishes  the  inertia  terms  to  be  used  based 
upon  the  keyed  input. 
S90    EO  =  EOi)  +  P  +  ^n 

910  F:>!=n:iij-p+2+(G+p*i^0) 
9;?o  t-:3=G+p^^o 
930   f:^(=xio-ki 

940    Fl=l+F-01 

9S0    F^=F-P*R0 

960    F0=1+F 

970     F3=G+F 

980    E5=E0*F1 -Fl  *H:'-^ 

990    E6=F0*F0-E1 *E1 

1000    F4=D1 *R1 +XIO+FO 

1010    F5  =  -^l  i<F!?  +  R*FO 

10^0    F6=E4*D1-XI0*F1 

1030    F7  =  D1-XI0*V/ 

1040    F8=F0--ei*V 

lObO    C  =  F1+F.4*V 

1060    C0=E3*(FH*E0-F1*E1 ) 

1070    C1=-EP*( (V)**. 5)*E3*F8 

10  80    C2=F7*(E0*C-E1 *F2 ) +E3*F8*(F2-R*V ) 

1090    C3=    -EP*(  (\y)**.  5)*F7*C 

1100    C4  =  F7*C*(F3-R*\/) 
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1110    C5=(F7*(E3  +  E1  )  )/C4 

1120    n6=(El*E3)/C4 

1130    G7=(F7*(F1 +E4+V) )/C4 

1140    C8=<E3*F8)/C4 

1150    C9=(E1 *F7+E3*F8)/C4 

1160    PRINT     100 

1170    100    Fgt^MATC*  C0,C5  C1*C6  02, C7 

1180    &C3*C8  C4,C9") 

1190    PRINT    1  10,C0,Cl,C2,C3,C4,C5,C6,n7,C8,C9 

1200    110    FO-^MATCIPSEM.S) 

1210    >J  =  4 

1220    CHK=0 

1230    Ad  )=C0;A(2)  =  0iA(3)=C2iAC4)=0;AC5)=C4 

1240    Aid )=0JAI(2)=C1 ;AI(3)=0iAI(4)=C3;AI(5)=0 

1250    120    CALL    CD0  WNH(  A.  A  I  ,N,  RR,C>^  ) 

1260    140    PRINT    150  ' 

1270    150    F0RMAT(/"  REAL    PART  IMAG.    PART 

1280    &  REAL    PART  IMAG.    PART") 

1290    PRINT    160,RRd  ),CRd  )*RR(2)*CR(2) 

1300    160    FORMATC    R00T    1     "  1  P2E1  4.  7,",    R00T    2    •'1P2E14.7) 

1310    PRINT    170,RR<3)*CRC3),RR(4),CR(4) 

1320    170    F0RMAT("    R00T    3    " 1 P2E1 4. 7,",    R00T    4    "1P2E14.7) 

1330    PRINT    180 

1340    180    F0RMATC/"    THE    REC0N.    C0EFF.    0F    THE    P0LYN.    ARE:") 

1350    PRINT    190 

1360    190    F0RMAT(/"    RE/IM--C0  CI  C2 

1370    &  C3  04") 

1380    PRINT    1  10,Ad  )*A(2),A(3),AC4),A(5) 
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Lines  890-1380  direct  the  calculation  of  the  coefficients,  the  solution  of  the 
Characteristic  Equation  for  the  roots  and  the  printing  of  the  reconstructed 
coefficients. 

1390    IF(EP)1  100,  n00#f?05 
1400    205    CONTINUE 

1410    PRINT     1 10#AI(1  ),AI(2),AI<3)*AI(4),AI(5) 
1420    1100     IFCKEY2-R)200*1000, 200 
1430    1000    G0T3(1O1O*1O2O* 1030*1040), L 

1440     1010    G3T0(79O,79O,79O,79O,46,44,42,44,44>,KEY1 
1450    1020    G0T0(79O#79O, 790, 790, 790, 790, 790,42, 42), KEYl 
1460    1030    G3T0(79O, 790, 790, 790, 790, 790, 790*790, 46), KEYl 
1470     1040    G0T0     790 
1480    200    055    204    K=  1  ,  4 
1490    RR\ (K)=RR(K) 
1500    204    CRl CK)=CRCK) 
1510    00    700    .J=l,51 
1520    9( J)=-1.02+.02*J 
1530    00    500    K=l,4 

The  logic  in  lines  1390-1530  uses  the  keyed  input  to  direct  control  to  the 
desired  location. 
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154OCCALCULAri0N  0F  ALPHA(K) 

1550  IF  (J.GT.1)G0T0  209 

1560  1=1 

1570  0^  206  -^=1,4 

1580  IF(K-M)208*206,208 

1590  208  ORR(  I  )  =  RRC'-<)-.^R(M) 

1600  OC^{CI)=C'^(K)-CK(M) 

1610  1=1+1 

1620  206  CONTINOE 

1630  AR  =  0'^-?(l  )  JXAI  =  0CK<1  )  ;BR  =  ORR  (2  )  ;8  I  =  DCR  C2  ) 

1640  XCR  =  ^ULT<AR#3R)-MULT(XAI*BI);C1=M1JLT(AR>BI)+MULT(XAI,-3R) 

16  50  AR=XCRiXAI=CIi3R=DRR(3>iBI=DCR(3) 

1660  XCR=MULT(AR,BR>-MULTCXAI,BI);CI=MULT(AR*BI)+MULT(XAI,BR) 

1670  AR=i .o;xai=o;br=xcr;bi=ci 

16  80  XCR=DI\/(AR*XAI,BR#BI);CI  =  0IV(XAI,-AR*BR>BI) 
1690  RACK)  =  XCRJCA(K)=C  I 

If  the  deflections  or  stresses  are  requested,  the  coefficient  a,   (used  in  the 
calculation  of  the  partial  fraction  expansion  of  the  Characteristic  Equation)  is 
determined  in  lines  1540-1690. 
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1700n    CHECK    T.3    SEE     IF    S0LIJTI0M    SATISFIES    THE    OIFF.     EQUATI3NS 

1710    Ar{  =  RR(K);XAI=CfR(K)JBR  =  RR(K)if3I=CRCK) 

1720    XCR=MULT(AR,BR)-MULT(XAI,BI)JCI=MULT(AR»BI)+MULT(XAI,^R) 

173  0    RMi|'v1=XC^lC^JUM=CIlAR=RA(K)  JXAI=CACK) 

17  40    XCR=DIV/(A*^,XAI,BR,9I);CI  =  0I\/(XAI*-AR,BR,RI) 

17  50    AR  =  E0+C7+C4  +  RNU'«1  +  C0J  XA  I  =  E0*C7*C4*CNUri 

1760    Bse=XCf^J9I  =  CI 

1770    XC-l  =  MlJLT(AR,RR)-M'JLT(XAI,8I>iCI  =  MiJLT(AR,BI  ) +-^iJLT  C  XA  I  ^'-S^  ) 

1780    CHK=XC^+CHK 

17  90     IF(K.NE. 4)G9T0    6001 

IROO    IJLIM  =  ABS(CHK+C4);  I  F  (  UL  IM-EPS  IL0N  )  600  1  *  6003  ,  6003 

IRIO    6003    PRINT    6004 

1820    6004    FP1RMAT(/"    THE    D-E.     IS    N3T    SATISFIED    BY    THE    SOL.") 

The  expression  derived  in  Appendix  G,  Equation  (G-11)  is  checked  in  the 
above  logic  to  ascertain  that  the  solution  satisfies  the  Differential  Equations. 
183OCCALCULATI0N    0F    BETACO 
1840    6001    C0'MTINIJE 

1850    AR=RR(K);XAI=CR(K);BR=RR(K)iBI=CR(K) 

1860    XCR  =  MULT(AR,BR)-MIJLTCXAI,BI);CI  =  MULT(AR»BI)+MULTCXAI,BR) 
1870    AR=(C7*XCR)+C8 
1880    XAI=C7*CI 
1890    BR=RA(K) JBI=CA(K) 

1900    XCR=iMULT(AR*BR)-MULT(XAI,BI)JCI  =  MULT(AR*BI)+MLJLT(XAI>BR) 
1910    R>3(K)  =  XCK 
1920    CB(K)=CI 
1930    209    CPItNjriNJlJE 

1940C    CALCULATION    0F    E    T0    THE     IS(K)0    AND     IS(K)C0+D)     P0WERS 
19S0    AR=    -0C J)*CR<K)iXAI=0C J)*RR(K) 
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1960  r)F:LrA  =  -30 

1 

1970  IF(Ai^.LT.nF"LrA)G0T'3    511 

1980  IFCAR.GF.  1  )G3  rP)21  1 

1  9  90  XCrJ  =  MIJLTCEXP(AK)>CGSCXAI  )  )  ;C  I  =  vllJLT  (EXP  (  AR  )  *  S  I^J  (  XA  I)  ) 

2000  GO  r-^    21? 

2010  211     Er>lK(K)  =  OJEPl  I  (K)  =  0 

2020  GO     ra     217 

2030  212    EP1R(K)  =  XC^^;EP1  ICK)=CI 

2040  217    C0iNiTIiMUE 

2050  AR=    -  (OC  J)+r))*CRCK)i     XA  1=  (0  <  J) +D  )  *:^R  (K  ) 

2060  IFCAR.Lr.DELTA)G0T0    213 

2070  IFCAR.GT. 1 )G0T0    213 

20  80  XCR=f'1ULTCEXP(AR),C0SCXAI)>;CI  =  VIULT(EXP(AR)>Sr^JCXAI)) 

30  90  G'3T0    214 

2100  213    EP2R(K)=0;EP2I <K)=0 

21 10  G0T0    218 

2120  214    EP2RCK)=XCR*EP2I(K)=CI 

2130  21 8    CONTINUE 


The  /3,   and  exponential  terms  are  calculated  above. 
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2140C    CftLCIILATI0N    OF    SG.M(BCK))     AND    THE    HEAVISIOE    FlJ,\)CTigMS 

■PA'^0  IFC.J.nr.  1  )G3  TP)     32S 

2160  A\y  =  A3SCCW(K)  ) 

2170  IFCAV.LT.  1  H:-4)CR<K)=0 

2180  295    IF(CR(K)  )300»310*320 

2190  310     IF(L3.FQ.O)G0    T(3    311 

2200  03    318    L4=l,4 

2210  RR(L4)='^R2(L4) 

2220  318    CR(L4)=CR2CL4) 

2230  0070    288 

2240  31 1    M=l 

2250  EP=SMALLV 

2260  G0T0    270 

2270  260    EP=10*Z 

2280  270    CI =-EP*( (V)**. 5)*E3*F8 

2290  C3=-EP*( (V)**. 5)*F7*C 

2300  AIC1)=0;AI(2)=C1;AI (3)=0;AI<4)=C3iAI(5)=0 

2310  CALL    CD3WNHCA,AI,N,RR»CR) 

2320  L3=l+L3 

2330  D0    273    L2=l*4 

2340  RR2(L2)=RRCL2) 

2350  273    CR2(L2)=CRCL2) 

2360  PRINT    275*M,EP,RR(1 ),CR(1 )#RR(2)»CRC2),RR<3),CR(3)# 

2365  ^.RRC4),CR(4) 

2370  275    F0RMATCI2*F7.4, 1P4E1 5. 7/9X# 1P4E1 5.7) 

2380  288     IF(CR(K) )300>290,320 

2390  290    M=M+1 iZ=EPiG0T026O 

2400  300    SGNCK)=-1 
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2A\0  G'lJ'l    325 

3420  320    SGN(K)=1 

2430  325     IFCSGN(K) )330,360,360 

2440  330     IF(T( J) )340,350>350 

2450  340     HO=l i     HgN=0 

2460  G3  rg    3  90 

2470  350     H3  =  0i     H^iM=l 

24  SO  GCITP)     3  90 

2490  360     IKCT (J) )370,380j3R0 

2500  370    H0=O;     H0M=1 

2510  G0Tg    3  90 

2520  380    H0=i;     H0iM  =  O 

2530  390     IF (SGN (K) ) 400, 430*430 

2540  400     IF(0CJ)+D)41O, 420*420 

2550  410    H3D=1;H0DM=O 

2560  G0T0    460 

2570  420    H0D  =  O;H0DN=1 

2580  6010    460 

2590  430     IF(0CJ)+D)44O*45O*45O 

2600  440    H0D=OJH0DN=1 

2610  G0TO     460 

2620  450    H0D=1  ;H0DiN  =  O 

2630  460    D0    465    Ll=l>4 

2640  RRCLl )=RR1 (LI ) 

2650  465    CRCLl )=CR1 <L1 > 

The  sgn  function  and  the  Heaviside  step  functions  are  determined  in  steps 
2140-2650. 
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2660C    CALCIILAIE    F(K)     AND    F(i<)PRIME 

2670    FPlr{CK)  =  M'^*EPlK<K);     FP  1  I  (K  )  =  HC5*EP1  I  (  K  ) 

2680    FP2'i:(K)  =  HgD*EP2K<K)     J  FP2  I  CK)  =  H0D*EP2  I  <K  ) 

2690  FPR(K)=FP1R(K)-FP2R(K) 

2700  8K=FPR(K) 

2710  FPI(K)=FP1 ICK)-FP2I(K) 

2720  9I=FPI(K) 

2730  FLH«CK)=SrrNJ  CK)*C9*CA(K) 

2740  AK=FLR(K) 

2750  FLI (K)=-SGN(K)*C9*RA(K) 

2760  XAI  =  FLI('r<) 

27  70    XC'^  =  MULTCAR*BR)-MULT(XAI,BniCI  =  '«lULT(ARjBI)+MIJLTCXAI,3R) 

2780  FR(K)=XCR;  FI(K)=CI 

27  90  AR=-CR(K);XAI=RR(K);3R=FR(K);BI=FI CK) 

2800  XCR  =  MULT(AR,RR)-MIJLT(XAI*8I)iCI  =  MULT(AR*BI)+MULT(XAIjBR) 

2810  nFR(K)=XCR;  DFI(K)=CI 

2820C  CALCULATIOM  f^F  WF(K) 

2830  WFPlR(K)=(-.5+EPlR(K) )*H0  + ( . 5 ) *H0N 

2840  WFPl I (K)=EP1 I (K)*H0 

2850  WFP2RCK)=C-. 5+EP2R(K))*H0D+(.5)*H0DM 

2860  WFP2I<K)=EP2I(K)*H0D 

2870  WFPRCK)=WFP1R(K)-WFP2R(K) 

2880  WFPI (K)=WFP1 I CK ) - WFP2 I (K) 

2890  AR=RB(K);XAI=C3(K) JBR=RRCK)iBI=CRCK) 

2900    XCR  =  DIV/CAR*XAI>BR*BniCI  =  DIV(XAI*-AR,BR*BI) 

2910    AR  =  SGN(K)*XCRJ  XA I  =  SGN (K > *C I ; BR  =  WFPR (K ) ; B 1  =  WFP I (K) 

2920     XCR  =  MULT(AR,BR)-MULT(XAI»BI  )iC  I=MIILT  (  AR,  B  I  ) +MUL  F  (  XA  I  *BR  ) 

2930    WFR(K)=XCRi     WFI(K)=CI 

2940    AR=R3(K);XAI=CB(K)1BR=FR(K)IBI=FI <K) 
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29b0    XCK  =  'yiULr(AR#BR)-i'^ULT(XAI*Bl)  JCI  =  ^^IJLT(A^,3I)+>^ULr(XAI,i;!'"0 

;?9A0    BK  =  -C9*RA(K)  ;  B  I  =  -C  9*CA  (K  )  i  AR=  XCR  S  XAI=CI 

3970    XCR  =  0IV(AR>XAI,BR*BI);CI  =  DW(XAI*-AK,B^»8I  ) 

29B0    DWFR(K)=XCRi     DWFI(K)=CI 

2990C    CALCULATION    0F    UF<K)    AND    THE    DERIVATIVE    0F    UFCK) 

3000    UFP1R(K)=(-.5*{1-CR(K)*0CJ))+EP1R(K) )*H3 

3010    &    +C.  5)*(1-CRCK)*0(.J)  )*H3N 

3020    UFPl I CK)=(-. 5*RR(K)*3CJ)+EP1 T (K)*HO+ 

3025    *     .5*(l+RK(K)+0(  J))*H0!M 

30  3  0    (iFP2R(K)=C-.  5*  (  1  -CR  (  K)  *C3  C  J) +0  )  ) +EP2R  (K  )  ) +HOD 

3040    &    +(.  5)*C1-CR(K)*(0(  J)+D)  )*H0DiN 

30  50    UFP2I(K)=(-.5*RR(K)*<0CJ)+D)+EP2I(K))*H0D 

3060    &    +.5*RRCK)*(0< J)+O)*H0DN 

3070    UFPR(K)=UFP1R(K)-UFP2RCK) 

3080    UFPI(K)=UFP1 I (K)-UFP2I (K) 

3090    BR=RR(K) JAR=RR(K)iBI=CR(K)iXAI=CR(K) 

3100    XCR=MULT(AR,BR)-MULT<XAI#BniCI  =  MULT(AR*Bn+MULT(XAI,BR) 

3110    AR  =  -SGN<K)*C6*CA(K);     XA  I  =  SGN (K ) *C6*RA (K ) ; BR= XCRiB I=C I 

3120    XCR=DIV(AR*XAI*BR*BI);CI=DIV(XAI,-AR,BR»BI) 

3130    UFLR(K)=XCRiUFLI (K )=C I ; AR=XCR J XA I=C I i8R=UFPR CK ) ; 

3135    &    BI=UFPICK) 

3140    XCR=MULT(AR,BR)-MULT(XAI»BI);CI=MULTCAR,BI)+MULT<XAI*BR) 

3150    iJFR(K)  =  XCRi     UFI(K)=CI 

3160    AR=RRCK);XAI=CRCK)iBR=RR(K)l8I=CR(K) 

3170    XCR=MULT(AR,3R)-MULT(XAI,BI)iCI=MULTCAR,BI)+MULTCXAI*BR) 

31 80    AR=WFRCK);XAI=WFICK);BR=C7*XCR+C8;BI=C7*CI 

31  90    XCR  =  DIV(AR#XAI*BR*BI);CI  =  Diy/(XAI*-AR*BR,BI) 

3200    DUFR(K)=-C6*XCRiDUFI(K)=-C6*CI 
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The  F,   functions  and  their  derivatives  are  calculated  in  lines  2660-3200. 
k 


:v->10  500    OONJTINUE 

32;?0C    CALCULAngN    0F    DISPLACEi'^lE.NJTS 

3230  PgORC  J)  =  F"RC1  )+FR<2)+FR(3)  +  F>^(4) 

3^40  P09ICJ)=FI(1 )+FIC2>+FI(3)+FI(4) 

32S0  W'30>^(.J)  =  l^JF,R(l  )+WFR(2)+WFR(3)  +  WFR(4) 

3260  W0CH  ( J)  =  WFI(1  )+WFI(2)+WFI <3)+WFI<4) 

3270  U'ZI^ORC  J)  =  -  (C5*P0QR(  J)  )/C9  +  UFK(  1  ) +UFR  (2  ) +  UFR  (  3  ) +  UFR  (  4  ) 

32  80  UOi)I(.J)  =  -  (C5*P0QIC  J))/n9  +  UFI  CI  ) +i.)FI  (2  ) +IJF  I  C3)+'JFT  (4) 
3290r;  CALCULATI0N  0F  FHF  DFRIVATIVES  0F  THE  0 1  SPLACE■^E^JT.S 
3300  DP  =  OFr^(  1  )+DFR(2)+DFRC3>+DFR(4) 

3310  DW    =r)WFU  1  )  +  0WFR(2)+DWFR(3)+DWFR(4) 

3320  niJF  =  DUFR(l  ) +DUFR  (2  )+DUFR  (3  ) +OUFR  (4  ) 

3330  DU=-(C5*DP)/C9+0UF 
3340C    CALCULATI0N    3F    THE    STRESSES 

33  50  SXMC J)=ElU*W3gR( J)+D1*DP+0U 
3360  STi'^C  J)  =  E0U*W3OR(  J)+E1U*DU 
3370  SXBC J)=.5*H*C0U+0P) 

33  80  STB(  J)  =  (.  5*H)*(E1  U*DP-ET0EX*W0G)R  (  J)  ) 

3390  SXZ(J)=G*CP0QR(J)+DW) 

3400  SX0( J)=SXM< J)-SXBC J) JST0 ( J)=STMCJ)-STB( J) 

3410  SXI(J)=SXM(J)+SX8(J)iSTI(J)=STMCJ)+ST9(J) 

3420  700  C0,MTINUE 

3430  G0T0(5O5*5O5, 505, 505*540, 540. 580)*KEY2 

3440  505  00  522  J=l,51 

3450  IFC J.HT. 1 )G0T0  515 

3460  PRINT  510 
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3470    510    FORMATC/    6X* 1 H0 > 1 5X, 1 HW,20X, 1 HU, 1 RX, 3HPS  I  ) 

3480    515    Pr^IiNJT    52O*3CJ)#W0QR(  J),IW9R(  J),P(?)OR(  J) 

3490    520    F0RMATC1F1O. 5, 1P3E20.6) 

3500    522    C!3^JTI^IUE 

3510    G9T0(1OOO,54O*54O*58O)*KEY2 

3520    540    DO     575    J=W51 

3530    IFC.J.GT.  1  )G0T3    560 

3540    PRIMT    550 

3550    550    Fg;^>lAT(/3X,  1  H0*8X,3HSXM»1OX*3HSXB,  10X*3HSTM,  lOX, 

3560    &3HST3> 10X>3HSXZ) 

3570    560    PRI^JT    570*  0  (  J)  *  SXM  (  J)  ,  SXB  (  J)  *  STMC  J)  *  STR  (  J)  ,  SXZ  C.J) 

3580    570    F0RMATC1F6.2,1P5E13.6) 

3590  575  C0NTINUE 

3600  G0T0(1OOO* 1000, 5 80* 5 80 #5 80, 1000),KEY2 

3610  580  00  610  J=l,51 

3620  IFCJ.GT. 1 )G0T0  600 

3630  PRIMT  590 

3640    590    F0RMAT(/3X, 1H0,8X,3HSX0,1OX,3HST0, 10X,3HSXI, lOX, 

3650    &3HSTI, 10X,3HSXZ) 

3660    600    PRINT    570,0 < J) , SX0 ( J), ST0 ( J) , SXI ( J) , STI ( J) , SXZ( J) 

3670    610    C0INJTINUE 

3680    G0T0     1000 

3690    790    PRINT    800 

3700    800    F0RMATC///////) 

3710    810    G0T0    5 

3720    820    ST0PJEMD 

The  displacements  and  stresses  are  calculated  and  printed  as  requested  and 
control  is  diverted  to  EMDATA  to  pick  up  other  cases  as  directed  in  lines 
3210-3720. 


APPENDIX  I 
RELATIONSHIPS  BETWEEN  ANALYSIS  PARAMETERS 

Dimensionless  and  Physical  Parameter  Relationships 
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For  external  pressure 

h       =     I  (1-9) 

O  2 

r       =     iH^  (I-IO) 


Isotropic  Material 

Hooke's  Law  for  a  homogeneous  isotropic  material  can  be  written 
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assuming  the  normal  stress  is  zero.    Comparing  these  expressions  with  the 

assumed  stress-strain  relationships  given  by  Equation  (2-10) 
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